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Executive  Summary: 

The  United  States  Navy  has  recently  embraced  the  revolutionary  eleetric  warship  coneept  for  the  design 
of  future  vessels.  A  core  element  of  this  eoneept  is  the  integrated  engineered  plant  (IEP).  1EP  removes 
traditional  system-level  barriers  between  the  distinet  ship  plants  (eleetrieal,  propulsion,  eooling, 
ventilation)  so  they  ean  leverage  a  shared  set  of  resources  and  information  management  systems.  The 
benefits  of  the  IEP  framework  are  naval  vessels  with  reduced  costs,  smaller  erews,  enhanced  effieieney 
and  improved  surv  ivability.  As  part  of  the  IEP  framew  ork,  plant  components  are  becoming  empowered 
with  computing  capabilities  for  self-control  and  management.  As  a  result,  the  IEP  is  a  complex  and 
highly  decentralized  system  that  is  still  responsible  for  achieving  real-time  control  objectives.  To  achieve 
optimal  control  within  the  challenging  environment  that  the  lEP’s  architecture  poses,  an  integrated 
hardware-software  framework  is  proposed.  At  the  core  of  the  proposal  is  the  design  of  a  modular  wireless 
node  prototype.  The  wireless  node  integrates  wireless  radios  and  embedded  computing  technologies  with 
existing  plant  sensors  and  actuators.  The  units,  termed  Narada,  are  low-eost  (costing  less  than  $200 
each),  eompaet,  hardened,  and  are  capable  of  eolleeting  sensor  data  and  issuing  actuation  commands. 

Wireless  sensor/aetuator  networks  forego  the  high  data  transfer  rates  associated  with  cabled  sensors  in 
exchange  for  low-eost  and  low  -power  communication  between  a  large  number  of  sensing  deviecs,  each  of 
whieh  features  embedded  data  processing  capabilities.  As  such,  a  newr  paradigm  in  large-scale  data 
processing  has  emerged;  one  where  communication  bandwidth  is  somewhat  limited  but  distributed  data 
processing  centers  are  abundant.  By  taking  advantage  of  this  grid  of  computational  resources,  data 
processing  tasks  once  performed  independently  by  a  eentral  processing  unit  ean  now  be  parallelized, 
automated,  and  carried  out  within  a  wireless  sensor  network.  By  utilizing  the  intelligent  organization  and 
self-healing  properties  of  many  wireless  networks,  an  extremely  scalable  multi-processor  computational 
framework  ean  be  developed  to  perform  advanced  engineering  analyses.  In  this  project,  a  novel 
parallelization  of  the  simulated  annealing  stochastic  search  algorithm  is  presented  and  used  to  update 
models  of  ship  plants  by  comparing  model  predictions  to  experimental  results.  The  resulting  distributed 
model  updating  algorithm  is  capable  of  being  used  as  a  computational  tool  to  perform  an  assessment  of 
the  condition  of  a  plant.  Should  damage  occur  to  a  plant,  it  ean  be  identified  in  a  timely  fashion  and 
reported  to  a  supervisory  element  in  the  plant  that  can  reconfigure  the  system  to  ensure  the  plant  ean 
continue  to  meet  its  mission  objective.  The  project  focused  on  the  monitoring  and  condition  assessment 
of  a  ehillcd  water  plant.  Using  a  demonstrator  located  at  the  Naval  Surface  Warefare  Center  (NSWC)  - 
Philadelphia,  the  Narada  network  proves  effective  in  identifying  rupture  conditions  in  the  plant. 

A  total  of  24  high  quality  publications  were  published  based  on  the  work  conducted  in  this  project.  In 
addition,  a  variety  of  technology  transfer  activities  were  undertaken  including  a  long-term  deployment  of 
the  wireless  nodes  developed  on  this  project  to  monitor  the  hull  of  the  FSF-1  SeaFighter,  a  high-speed 
aluminum  littoral  combat  ship.  During  summer  2009,  a  31 -channel  wireless  sensor  network  w'as  installed 
to  reeord  hull  strains  and  accelerations  with  the  ship  monitored  as  it  travelled  from  Panama  City  to 
Portland,  Oregon. 
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1.  Introduction: 

To  maintain  its  role  as  the  preeminent  global  naval  power,  the  United  States  Navy  invests  heavily  in 
fundamental  research  aimed  at  modernizing  its  fleet.  As  part  of  these  efforts,  the  Navy  has  begun  to 
explore  revolutionary  ship  design  concepts  intended  to  enhance  mission  functionality,  improve  damage 
tolerance,  and  reduce  manning  requirements.  In  recent  years,  the  electric  warship  design  paradigm  has 
emerged  as  one  such  revolutionary  concept  (ONR  2003).  The  electronic  Inaval  warship  seeks  to  employ 
excess  power  generated  in  the  shipboard  propulsion  system  to  power  auxiliary  ship  systems  including 
power-pulsed  weapons,  high  energy  embedded  sensors,  among  other  systems.  Benefits  associated  with 
Integrated  Power  Systems  (IPS)  include  reduced  cost,  enhanced  naval  design  flexibility,  and  increased 
mission  efficiency.  In  tandem  with  IPS,  the  Navy  has  also  pursued  enhancement  (namely  improved  ship 
reliability  and  survivability)  of  future  naval  vessels  through  the  Integrated  Engineering  Plant  (IEP) 
architecture. 

1.1.  Integrated  Engineering  Plant  Concept: 

The  IEP  architecture  seeks  further  integration  of  all  ship  plants  (e.g.,  propulsion,  cooling,  ventilation,  fire 
suppression,  electrical)  so  that  they  share  common  resources  and  information  management  systems 
(Dunnington,  et  al  2003).  IEP  capitalizes  on  the  continuing  distribution  of  computing  elements 
throughout  the  ship.  As  microprocessor  electronics  become  increasingly  capable  (from  a  computing 
standpoint)  and  inexpensive,  computing  resources  are  moving  away  from  centralized  computing  systems 
(as  was  typical  of  older  naval  plant  designs)  and  moving  towards  integration  w  ith  the  sensor  and  actuators 
embedded  within  the  distinct  shipboard  plants.  The  trend  of  integrating  computing  resources  with 
component-level  sensor  systems  allows  the  IEP  architecture  to  adopt  automated  decision  making  in  a 
decentralized  fashion. 

At  the  core  of  IEP  is  a  three  tiered  decision  and  control  architecture  divided  between  the  component, 
process  and  mission  control  layers,  as  shown  in  Figure  1.1  (Dunnington,  et  al.  2003).  The  component 
layer  consists  of  the  various  machines  and  modules  in  which  sensors  and  actuators  reside.  Some 
computing  technology  resides  within  the  component  layer  to  carry  out  tasks  associated  w  ith  local  closed- 
loop  control  and  component-level  health  monitoring.  The  process  layer  resides  above  the  component 
layer  and  contains  greater  computing  resources.  Given  the  responsibility  of  optimizing  the  performance 
of  the  different  ship  objectives,  the  process  layer  continuously  monitors  the  component  layer  to  be  aware 
of  the  need  to  reconfigure  ship’s  resources  to  better  attain  the  ship’s  mission.  The  highest  layer  of  the 
three  tier  hierarchy  is  the  mission  control  layer.  The  mission  control  layer  manages  the  global 
performance  of  the  ship  by  making  major  resource  allocations  when  the  operating  mode  and  physical 
condition  of  the  ship  changes.  At  this  layer,  the  IEP  system  provides  ship  personnel  with  critical  real¬ 
time  information  collected  from  the  lower  component  and  process  layers. 


Figure  1.1:  Integrated  engineering  plant  (IEP)  concept. 
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The  mission  efficiency  of  the  ship  is  drastically  improved  because  of  the  shared  and  interoperable 
management  system  provided  by  the  IEP  framework.  For  example,  as  sensor  and  plant  systems  continue 
grow  in  complexity  and  become  increasingly  interdependent,  automated  decision  making  can  alleviate  the 
problem  of  inundating  the  ship  crew  with  too  much  information  to  be  effective  in  their  respective 
missions.  Hence,  automated  decision  making  can  lead  to:  1)  substantial  reductions  in  the  size  of  ship 
crews,  2)  improvement  of  the  mission  responsiveness  of  the  ship,  and  3)  rapid  reallocation  of  plant 
resourced  for  ship  survivability.  A  highly  decentralized  system  also  derives  benefits  by  avoidance  of 
single-point-of-failure  problems  common  in  centralized  control  system  architectures.  Some  of  the  other 
benefits  offered  by  the  IEP  framework  include  system  scalability,  cost-effcctiveness,  health  monitoring  of 
shipboard  equipment  and  increased  payload  fractions. 

To  render  the  IEP  framework  feasible,  a  reliable  ship-wide  area  network  (SWAN)  is  needed  to  provide  a 
communication  infrastructure  within  and  between  the  three  IEP  layers.  The  required  SWAN  must  be 
immune  to  physical  damage  that  is  possible  during  battle.  Current  ship  designs  eall  for  triple  redundant 
serial  buses  for  communication  between  components.  In  contrast  to  this  approach,  various  military 
vendors  have  proposed  different  communication  architectures  that  enjoy  better  fault  tolerance.  For 
example.  Adept  Systems  has  proposed  a  complex  hybrid  communication  architecture  comprising  of 
multi-point  ring  networks  and  a  LonTalk  communication  protocol  (ANSI  709.1)  (Smith  1994).  This 
approach  is  scalable,  cost-effective  and  enjoys  self-healing  capabilities  as  recently  illustrated  on  the  ONR 
Afloat  Laboratory  test-bed  ship  and  the  cx-USS  Peterson  (Adept  2004).  Highly  reliable  network 
architectures  are  also  being  explored  for  the  distributed  power  systems  of  future  naval  vessels.  For 
example,  ring  topology  load  buses  sharing  strong  similarity  to  the  multi-point  ring  communication 
networks  proposed  by  Adept,  are  proposed  as  part  of  the  reconfigurablc  electrical  systems  (RES)  program 
(Dunnington,  et  al.  2003). 

1.2.  Wireless  Communications  in  Naval  Ships: 

As  the  United  States  Navy  progresses  towards  the  Integrated  Engineering  Plant  design  concept,  the 
extensive  length  of  wires  that  are  needed  to  connect  sensors  and  actuators  with  control  units  distributed 
throughout  a  vessel  remains  a  technological  challenge.  In  addition  to  the  high  cost  of  routing  wires 
during  construction,  naval  vessels  represent  a  complex  and  harsh  environment  in  which  extensive  lengths 
of  wires  arc  vulnerable  to  detriments  such  as  heat,  moisture  and  toxic  agents  (MacGillivray  and  Goddard 
1997).  With  wires  vulnerable  to  failure  when  exposed  to  these  harsh  conditions,  reduction  or  outright 
elimination  of  communication  wires  would  greatly  enhance  the  reliability  of  on-board  engineering  control 
systems  in  addition  to  reducing  installation  and  maintenance  costs. 

Ships  pose  a  challenging  setting  for  the  propagation  of  wireless  signals.  Estes,  et  al.  (2001)  has  explored 
the  feasibility  of  wireless  radios  for  both  intra-  and  inter-eompartment  shipboard  communications  within 
various  naval  vessels  (cx-USS  America,  USS  Ross,  and  USS  Carr).  Their  study  considered  radio 
frequencies  between  800  MHz  and  2.5  GHz  which  arc  typical  radio  frequencies  for  eommereial  off-the- 
shelf  (COTS)  radios.  Based  on  the  high  radio  frequency  (RF)  reflectivity  of  steel,  multi-path  effects  were 
discovered  to  dominate  received  radio  signals  during  inter-compartmcnt  wireless  communication.  To 
overcome  multi-path  influences,  only  frequency  hopping  spread  spectrum  (FHSS)  wireless  radios  were 
discovered  to  effectively  work.  When  assessing  the  feasibility  of  inter-compartmcnt  communication, 
ship  bulkheads  severely  attenuated  wireless  signals  (on  the  order  of  magnitude  of  20-30  dB)  but 
communication  through  two  or  three  bulkheads  was  still  found  to  be  possible.  Steel  is  a  near  perfect 
conductor  that  reflects  electromagnetic  waves  thereby  limiting  radio  signal  penetration.  However,  on 
modem  ships,  a  number  of  non-steel  elements  arc  present  in  the  bulkheads  (c.g.  hatch  seals,  ducts,  cable 
transits)  that  allow  wireless  signals  to  penetrate.  Mokolc,  et  al.  (2000)  has  undertaken  a  similar  wireless 
communication  feasibility  study  on  the  cx-USS  Shadwell  using  COTS  wireless  modems  that 
communicate  on  the  800  MHz  to  3  GHz  radio  frequencies.  Their  study  found  that  radio  communication 
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is  possible  for  inter-eompartmcnt  communication  using  commercial  wireless  radios,  even  when  bulkhead 
closures  were  fastened. 

In  1999,  the  Office  of  Naval  Research  (ONR)  initiated  the  Reduced  Ship’s  Crew  by  Virtual  Presence 
(RSVP)  program  with  the  aim  of  developing  a  proof-of-concept  wireless  monitoring  system  for  shipboard 
sensor  systems  (Schwartz  2002;  Seman,  et  al.  2003).  Wireless  communication  between  embedded 
sensors  and  existing  shipboard  local  area  networks  (LANs)  was  aimed  to  reduce  ship  construction  and 
maintenance  costs  in  addition  to  enhancing  the  reliability  and  surv  ivability  of  ships.  The  complete  RSVP 
framework  proposed  dedicated  sensor  networks  for  monitoring  the  environmental  parameters  of  ship 
spaces,  the  structural  integrity  of  the  hulk  the  health  of  machinery,  and  the  real-time  location  of  ship 
personnel.  The  framework  sought  to  minimize  the  power  consumption  of  sensing  devices  to  ensure  long 
operational  lives.  The  machinery,  environmental,  structural  and  personnel  wireless  sensor  networks  had 
redundant  access  to  various  access  points  through  802.11  wireless  communications.  Data  from  the 
distributed  wireless  sensors  are  communicated  to  a  dedicated  watch  station  through  a  highly  reliable  wire- 
based  shipboard  LAN.  The  RSVP  program  validated  different  components  of  the  RSVP  architecture  on 
numerous  naval  vessels  including  the  USS  Monterey  and  ex-USS  Shadwell. 

Other  researchers  have  successfully  explored  the  use  of  wireless  communications  integrated  within 
various  shipboard  engineering  plant  systems  to  monitor  component  health.  Ploeger,  et  al  (2003) 
describes  a  cost-effective  wireless  monitoring  system  that  monitors  the  operational  health  of  a  shipboard 
ventilation  system.  The  wireless  system  proposed  is  constructed  from  multiple  wireless  data  acquisition 
nodes,  termed  the  Intelligent  Component  Health  Monitor  (ICHM),  that  are  capable  of  collecting  sensor 
data  from  numerous  analog  sensors  and  communicating  that  data  via  Bluetooth  wireless  radios  to  a 
centralized  data  repository  termed  the  Compartment  Health  Monitor  (CHM).  The  system  described  is 
well  suited  for  intra-compartment  communication  because  of  the  short  20-30  m  communication  range 
associated  with  Bluetooth  radios.  Simple  data  processing  of  sensor  data,  including  threshold  detections 
and  spectrum  analysis,  is  performed  at  the  compartment’s  CHM  server.  As  validation,  ventilation  fans 
upon  an  operational  aircraft  carrier  were  monitored  for  overall  health  using  the  described  system, 
resulting  in  an  estimated  shipboard  crew  labor  reduction  of  63%  (reduction  of  5000-8000  man  hours  per 
year  per  aircraft  carrier)  (Ploeger,  et  al.  2003). 

Many  other  engineering  domains  have  begun  to  explore  the  use  of  wireless  Sensors  for  reducing  the  high 
costs  and  technical  difficulties  often  associated  with  wire-based  data  acquisition  systems.  For  the  past 
decade,  academic  and  commercial  research  teams  have  proposed  a  variety  of  low-cost  wireless  sensor 
platforms.  The  most  notable  examples  include  the  Smart  Dust  (University  of  Califomia-Berkeicy)  and 
pAMPS  (MIT)  projects  that  have  yielded  low-eost  wireless  sensor  nodes  designed  for  deployment  in 
wireless  sensor  networks  defined  by  high  nodal  densities  (Horton,  et  al  2002;  Min,  et  al  2001). 
However,  these  generic  wireless  sensing  systems,  in  their  current  state  of  development,  do  not  address  the 
unique  demands  of  the  harsh  monitoring  domain  posed  by  a  complex  naval  vessel,  where  low-power 
consumption,  far-reaching  communication  ranges,  penetration  characteristics  through  ship  materials  and 
sufficient  computational  capabilities  for  autonomous  data  processing  are  all  needed. 

1.3.  Decentralized  Large-Scale  Control  and  Artificial  Intelligent  Agents: 

As  industrial  systems  grow  in  size  and  complexity,  the  ability  to  control  them  in  centralized  fashion 
becomes  increasingly  difficult.  In  particular,  a  number  of  factors  contribute  to  the  breakdown  of  the 
centralized  control  system  architecture.  First,  as  the  number  of  sensors  and  actuators  linearly  increase,  the 
computational  demands  (aggregation  of  sensor  data,  calculation  of  control  forces,  and  commanding 
actuators)  placed  upon  a  centralized  controller  grow  exponentially  (Lunzc  1992).  Furthermore,  large 
plants  are  often  complex  systems  that  are  difficult  to  precisely  model  using  classical  mathematical 
approaches.  In  response  to  these  limitations,  the  field  of  control  engineering  has  extensively  explored  the 
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decentralization  of  the  control  system.  By  dividing  a  complex  global  system  into  simpler  sub-systems, 
optimal  closed-loop  control  can  still  be  achieved  by  decentralized  control  (Siljak  1991). 

Other  researchers  have  pursued  the  use  of  agents  for  performing  control  of  large-scale  complex  systems. 
Agents  are  software  abstractions  that  are  capable  of  perceiving  their  environment  (through  sensing)  and 
acting  upon  that  perception  (through  actuation)  based  on  a  rational  decision  process  (through  embedded 
intelligence)  (Jennings,  et  al  1997).  As  part  of  the  rational  decision  process,  the  agent  seeks  to  maximize 
a  performance  measure  that  represents  the  return  on  investment  for  a  given  action.  Agents  have  found  use 
in  a  number  of  complex  systems  including  control  systems,  networks,  software  systems,  manufacturing 
processes,  among  many  others  (Wooldridge  and  Jennings  1995).  Autonomous  intelligent  agents,  when 
assembled  into  a  community,  are  often  referred  to  as  multi-agent  systems  (MAS).  MAS  are  particularly 
attractive  for  performing  control  of  complex  systems  because  they  represent  a  scalable  problem  solving 
approach  for  problems  with  incomplete  information  and  decentralized  data.  For  control  problems  where  a 
centralized  control  solution  is  not  feasible  (due  to  the  inability  to  exactly  represent  the  system  with 
mathematical  models),  multi-agent  systems  have  proven  to  be  reasonable  approaches  to  performing 
feedback  control. 

The  1EP  framework  proposed  for  future  naval  ships  presents  many  technical  challenges  that  renders 
centralized  control  an  impractical  approach.  Alternatively,  the  use  of  multi-agent  systems  have  been 
proposed  as  a  means  of  controlling  the  1EP  in  a  more  effective  and  reliable  manner.  With  high  degrees  of 
autonomy  and  the  ability  to  rapidly  adapt,  intelligent  agents  will  undoubtedly  play  a  major  role  in  the  1EP 
framework.  Agents  provide  a  reconfigurablc  framework  for  plant  control  that  can  adapt  with  changes  in 
the  ship  mission  or  ship  operational  conditions  (an  undamaged  versus  damaged  ship).  The  agent 
framework  proposed  for  the  1EP  is  possible  largely  because  computing  power  is  migrating  from  upper 
tiers  of  the  IEP  framework  down  to  the  component  level.  Already,  researchers  are  validating  the 
feasibility  of  MAS  for  performing  adaptive  control  within  the  IEP.  Sun  and  Cartes  (2004)  propose  the 
use  of  software  agents  for  performing  closed-loop  control  shipboard  distributed  power  plant.  Maturana, 
et  al.  (2004)  describes  the  implementation  of  intelligent  agents  for  controlling  the  flow  of  water  within  a 
chilled  water  system.  Their  study  has  implemented  the  agents  using  legacy  control  modules  with 
abstracted  Java-based  software  agents;  validation  of  their  multi-agent  system  is  performed  upon  a  small- 
scale  chilled  water  plant. 


2.  Technical  Objectives: 

This  project  explores  the  use  of  wireless  sensor  and  actuator  networks  within  naval  vessels  as  a  means  of 
reducing  crew  requirements  and  to  enhance  ship  fight-through  capabilities.  Wireless  sensors  and  wireless 
actuators  are  not  sensors  or  actuators  per  se ,  but  rather  autonomous  nodes  of  a  data  acquisition  and 
feedback  control  to  which  sensors  and  actuators  could  be  attached.  The  advantages  of  using  wireless 
sensors  for  this  purpose  arc  to:  1)  reduce  the  installation  and  maintenance  costs  associated  with  current 
tethered  control  systems,  2)  to  allow  for  ad-hoc  network  topologies  between  sensors  and  actuators 
offering  communication  robustness  during  battle,  and  3)  to  collocate  computing  intelligence  with  plant 
components  for  plant  assessment,  control  and  reconfiguration.  This  represents  a  highly  revolutionary 
approach  to  the  implementation  of  the  integrated  engineering  plant  framework. 

2.1.  Objective  1:  Low-Cost,  Computationally-Rich  Wireless  Node  Development: 

The  first  technical  objective  of  the  project  is  the  design  of  a  low-cost  wireless  node  designed  explicitly  for 
embedment  in  shipboard  plant  systems.  Unlike  generic  commercial  wireless  sensors,  the  wireless  node  to 
be  developed  as  part  of  this  project  will  be  designed  to  achieve  functionality  associated  with  the 
monitoring,  control  and  reconfiguration  of  shipboard  plants.  Towards  this  end,  the  node  will  be  designed 
to  perform  the  following  tasks:  1)  collection  of  state  data  from  sensors  monitoring  plant  performance;  2) 
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real-time  wireless  communication  of  data  to  the  upper  ship  mission  tiers;  3)  embedded  data  interrogation 
for  intelligent  in-network  data  processing;  4)  ad-hoe  communication  between  wireless  sensors  to  transfer 
real-time  state  data;  5)  collective  decision  making  between  wireless  sensor  nodes;  5)  automated  actuation 
of  plant  elements  at  the  component  tier.  In  addition  to  these  features,  the  w  ireless  sensors  must  be  low- 
eost  (<$200)  and  be  hardened  to  ensure  they  survive  during  battle  and  when  exposed  to  water. 

2.2.  Objective  2:  Embedded  Firmware  for  Wireless  Node  Autonomy: 

Onee  the  wireless  sensors  are  designed  and  constructed,  the  sceond  project  objective  if  foeused  upon  their 
autonomous  operation.  To  automate  their  operation,  embedded  software  is  needed  for  eaeh  wireless 
sensor.  Embedded  software,  termed  firmware,  will  automate  the  collection  of  state  data  from  sensors, 
w  irelessly  eommunieate  data  with  other  wireless  sensors  and  to  issue  command  signals  to  plant  aetuators. 
A  eorc  element  of  the  firmware  is  the  eoding  of  a  flexible  wireless  communication  protocol.  Provided  the 
low-power  nature  of  the  1EEE802.15.4  communication  standard,  the  1EEE802.15.4  physical  and  medium 
aeeess  eontrol  protocol  layers  will  be  written  in  C  and  eneoded  in  the  wireless  sensor  eores.  With  a 
flexible  communication  standard,  the  wireless  sensor  network  will  be  able  to  form  ad-hoe  network 
topologies  between  plant  components. 

2.3.  Objective  3:  Distributed  Computing  for  State  Awareness  (Damage  Detection): 

The  third  objective  of  the  project  is  to  take  full  advantage  of  the  embedded  intelligence  eoupled  with  each 
wireless  sensor  through  formulation  of  agent-behavioral  models.  Borrowing  eoneepts  from  the  field  of 
artificial  intelligence  (AI),  eaeh  wireless  sensor  is  modeled  as  an  autonomous  agent  eapable  of  interaction 
with  other  wireless  sensors.  The  agent  framework  is  especially  powerful  for  achieving  decentralized 
eontrol  in  a  network  of  wireless  sensors.  Computationally,  the  project  foeus  will  be  on  the  development 
of  computational  tools  that  ean  perform  model  updating.  Using  an  analytical  model  that  describes  the 
configuration  and  operation  of  a  plant  system,  model  updating  ean  be  performed  to  realize  the  best  fit  of 
the  analytical  model  to  experimental  data  eolleets.  If  a  plant  experience  damage,  this  state  will  be 
identified  by  the  model  updating  method.  With  model  updating  inherently  an  intraetable  combinatorial 
optimization  problem,  stoehastie  methods  that  offer  high  probabilities  of  convergence  to  global  minima 
w  ill  be  adopted.  Furthermore,  to  benefit  from  a  dense  netw  ork  of  w  ireless  sensors  and  aetuators  installed 
in  a  ship  plant,  a  means  of  elegantly  distributing  computational  tasks  to  available  wireless  nodes  will  be 
developed.  Using  a  multi-agent  approach  based  on  free  market  economies  (termed  market  based  eontrol) 
will  be  explored  for  solving  eomplex  computational  resource  allocation  problems  encountered  in 
distributed  in-network  data  processing. 

2.4.  Objective  4:  Experimental  Validation: 

The  final  objective  is  to  perform  various  validation  experiments  of  the  integrated  hardware/software 
framework  proposed.  Opportunities  will  be  sought  to  validate  the  wirelessly  integrated  plant  on  testbed 
platforms  that  exist  at  NSWC  Carderoek  Division  -  Philadelphia.  Specifically,  a  scaled  down  chilled 
wratcr  plant  demonstrator  will  be  utilized. 


3.  Design  of  a  Low-Cost  Wireless  Sensor/Actuator  Node  for  Ship  Automation: 

3.1.  Design  of  a  Wireless  Node: 

While  wireless  sensor  solutions  are  commercially  available,  these  solutions  are  too  generie  for  use  in 
naval  plants.  A  more  effective  approach  is  to  design  a  wireless  node  whose  performance  has  been 
optimized  for  operation  in  a  naval  vessel.  A  wireless  sensor  will  be  designed  for  real-time  eontrol  of 
naval  plants  using  eommereial  off-the-shelf  (COTS)  components.  The  wireless  sensor  design  will  be 
broken  down  to  four  functional  subsystems:  sensing  interface,  actuation  interface,  computational  core  and 
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wireless  communication  ehannel.  The  sensing  interface  is  responsible  for  data  collection  from  plant 
sensors,  the  computational  eore  is  intended  to  aggregate  data  and  locally  process  it,  the  actuation  interface 
will  be  used  to  command  plant  actuators,  and  the  wireless  ehannel  will  establish  connectivity  between 
sensors,  actuators  and  the  ship  information  management  system.  An  architectural  schematic  of  the 
proposed  wireless  node,  termed  Narada,  is  shown  in  Figure  3.1  and  Figure  3.2. 

5.7.7.  Sensing  Interface: 

The  sensing  interface  is  the  wireless  node  component  to  which  plant  sensors  are  attached.  The  wireless 
node  should  be  able  to  aeeommodatc  multiple  sensors  at  one  time.  Furthermore,  the  interface  should  be 
independent  of  the  sensor  type  to  ensure  versatility  in  deployment  in  naval  systems.  The  primary 
eleetrieal  component  selected  for  the  sensing  interface  is  a  four-ehannel,  16-bit,  100  kHz  analog-to-digital 
converter  (ADC);  in  particular,  the  Texas  Instruments  ADS8341  is  selected.  The  interface  is  eapable  of 
accepting  analog  sensor  signals  that  are  between  0  and  5  V. 

3.7.2.  Computational  Core: 

The  wireless  node  has  at  its  eore  a  low  voltage  microcontroller  which  is  used  to  operate  all  of  the 
eleetrieal  components,  in  addition  to  performing  data  interrogation  tasks.  The  microcontroller  selected  is 
the  8-bit  Atmel  ATmegal28  microcontroller.  The  Atmegal28  provides  128  kB  of  flash  memory  which  is 
adequate  for  most  embedded  software  algorithms  required  for  monitoring  and  eontrol  applications.  To 
augment  the  built-in  4kB  of  SRAM,  external  memory  of  128kB  is  provided  using  a  Cypress  CY62128B 
SRAM  chip.  This  is  adequate  for  data  storage  and  analysis  requirements  of  the  sensing  unit.  The 
ATmcgal28  is  operated  at  8  MHz  whieh  consumes  less  than  100  mW  of  power  from  portable  power 
sources  ( e.g .,  batteries).  Power  is  a  consideration  because  it  is  anticipated  the  wireless  nodes  might  be 
required  to  run  on  batteries  during  some  battle  scenarios. 

3.7.3.  Wireless  Communication  Channel: 

A  low-power,  single-chip,  IEEE  802.15.4  compliant  wireless  radio  is  adopted  in  the  wireless  node  design. 
The  newly  created  IEEE  802.15.4  standard  is  especially  suitable  for  distributed  computation,  as  it  was 
developed  for  true  ad-hoe  peer-to-peer  networking  among  battery  powered  wireless  devices.  The 
Chipeon  CC2420  wireless  transceiver  is  selected.  The  CC2420  is  a  single  ehip  IEEE  802.15.4  compliant 
radio  eapable  of  ranges  adequate  for  intra-eompartment  communication  in  ships.  For  example,  its 
operational  communication  range  is  validated  up  to  20  m.  Furthermore,  the  radio  operates  on  the  2.4  GHz 
radio  band  with  data  rates  as  high  as  250kpbs.  Beeause  the  CC2420  is  designed  specifically  to  be  IEEE 
802.15.4  compliant,  it  does  not  require  the  use  of  a  central  coordinator  to  facilitate  communication;  the 
sensors  may  communicate  directly  with  each  other  making  them  ideal  for  peer-to-peer  communications. 

3. 1.4.  Actuation  Interface: 

The  Texas  Instruments  DAC7612  has  been  integrated  in  the  design  of  the  wireless  node  to  serviee 
actuators.  The  DAC  allows  two  actuation  channels  to  operate  simultaneously.  The  output  voltage  range 
of  the  DAC  is  0  to  4  V. 

3.2.  Assessment  of  Data  Collection  Quality  : 

With  a  set  of  wireless  nodes  fully  assembled,  various  diagnostic  tests  were  conducted  to  validate 
assumptions  made  during  the  design  proeess  of  the  node.  Specifically,  the  ADC  resolution  was  measured 
using  a  eonstant-voltage  battery  souree.  The  sensing  interface  resolution  is  an  important  property  that 
describes  the  data  eolleetion  aeeuraey  of  the  wireless  sensor.  For  example,  a  high-resolution  ADC  (16- 
bits  or  higher)  suggests  the  wireless  sensor  is  eapable  of  reading  sensors  with  low  amplitude  outputs.  In 
contrast,  a  low-resolution  ADC  (less  than  16-bits)  represents  a  eruder  data  eolleetion  tool  where  the 
quantization  proeess  of  the  ADC  ignores  small  varying  sensor  outputs.  Particularly  for  naval  plant 
systems,  sensors  defined  by  high  aeeuraey  must  be  complemented  by  a  high  resolution  ADC. 
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Figure  3.1  Design  of  a  wireless  node  for  ship  systems:  (a)  hardware  design;  (b)  nine  fully 
assembled  wireless  sensors  (a  blank  printed  circuit  board  is  shown  in  upper  left  corner). 


Figure  3.2.  Final  fully  assembled  wireless  sensor  (Chipcon  transceiver  not  shown) 


While  the  ADC  (Texas  Instruments  ADS8341)  has  a  nominal  resolution  of  16-bits,  the  presence  of  digital 
circuit  elements  ( e.g .,  microcontroller,  radio)  on  the  same  circuit  board  can  result  in  a  reduction  of  the 
ADC  resolution.  Digital  circuit  elements  internally  contain  electrical  gates  that  open  and  close  during 
their  operation.  As  these  electrical  gates  open,  electrical  charge  floods  into  the  wireless  sensor  ground 
plane.  If  the  analog  and  digital  circuits  share  a  common  ground  plane,  excess  charge  Hooding  can  result 
in  the  ground  exhibiting  an  apparent  change  in  voltage  relative  to  the  power  source  of  the  system.  The 
ADC  employs  two  references  that  are  assumed  to  be  fixed  to  convert  sensor  outputs  to  digital  formats:  a 
high  voltage  reference  (5  V  power  plane)  and  a  low  voltage  reference  (ground  plane).  Should  those 
references  change  in  any  way,  the  ratiometric  conversion  of  the  ADC  analog  input  (i.e.,  sensor  output) 
would  exhibit  quantization  noise  (random  toggling  of  the  least  significant  bits  of  the  ADC  digital  output). 
The  toggling  of  more  than  the  16th  least  significant  bit  of  a  16-bit  ADC  signifies  a  reduction  of  the  true 
ADC  resolution.  The  primary  motivation  of  pursuing  a  4-layer  circuit  board  design  w^as  to  employ 
separate  power  and  ground  planes  for  analog  and  digital  circuit  elements.  Even  with  a  4-layer  board, 
some  resolution  reduction  is  anticipated. 
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Figure  3.3.  Integer  output  of  each  wireless  sensor  channel  when  a  1.5  V  battery  is  attached. 

The  true  resolution  of  the  w  ireless  sensor  prototype  was  measured.  An  alkaline  battery  (1.5  V)  was 
interfaced  to  each  sensor  channel  of  the  wireless  sensor  (4  channels  in  total).  Batteries  are  ideal  power 
sources  for  such  tests  because  they  exhibit  noise- free  constant  voltages.  The  w  ireless  sensor  recorded  the 
battery  analog  voltage  as  a  16-bit  digital  number;  the  number  was  also  w  irelessly  transmitted  to  a  remote 
data  repository  w  here  time-histories  of  the  sensor  channel  voltages  were  stored.  The  time-history  plots  of 
the  measured  battery  voltage  are  show  n  in  Figure  3.3  for  each  of  the  4  sensor  channels.  The  time -history 
plots  are  in  terms  of  the  ADC  integer  values;  6000  data  points  collected  at  200  Hz  on  each  sensor  channel. 


To  quantify  the  true  resolution  of  the  ADC,  the  time-history  plots  of  the  ADC  integer  outputs  were  used 
to  generate  histograms  of  the  ADC  integer  outputs.  To  enhance  the  clarity  of  the  histograms,  the  mean 
integer  value  is  subtracted  from  each  time  history  record.  As  shown  in  Figure  3.4,  the  zero-mean  integer 
histograms  for  each  sensor  channel  exhibited  Gaussian  distributions  as  i$  typical  of  ADC  quantization 
error  processes.  The  majority  (>98%)  of  the  ADC  integer  values  fell  w  ithin  -4  to  4.  This  means  the 
integer  outputs  vary  by  8  integer  values  which  is  equivalent  to  3-bits  (8  *=  23).  If  the  resolution  of  the 
ADC  is  16-bits,  then  only  the  last  least  significant  bit  is  anticipated  to  be  metastable  (toggling  between  0 
and  1).  In  the  ease  of  the  wireless  sensor,  the  last  three  bits  w  ere  toggling.  As  such,  the  w  ireless  sensor’s 
true  ADC  resolution  was  assessed  to  be  14-bits. 
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Figure  3.4.  Performance  of  the  on-board  analog-to-digital  converter  (ADC)  -  histogram  of  ADC 
integer  output  to  assess  true  conversion  resolution. 


Figure  3.5.  (a)  Wireless  sensor  prototype  commanded  to  output  a  high-frequency  sinusoidal 
wave  on  its  actuation  interface  (wireless  sensor,  7.5V  battery  source,  and  Agilent  oscilloscope 
shown  in  this  laboratory  demonstration);  (b)  sample  actuation  signals  generated  by  the  wireless 
sensor  and  recorded  by  the  oscilloscope. 


3.3.  Illustration  of  Actuation  Interface  Operation: 

The  wireless  sensor  is  capable  of  outputting  0  to  4  V  analog  voltage  signals  using  its  on-board  2-ehannel 
DAC.  To  illustrate  the  operation  of  the  actuation  interface,  a  standard  laboratory  oseilloseope  is 
interfaced  to  the  actuation  interface  to  measure  the  output  signal  generated  by  the  wireless  sensor.  Figure 
3.5  presents  the  laboratory  set-up  of  the  wireless  sensor  interfaced  to  the  oseilloseope.  The  wireless 
sensor  is  commanded  to  output  a  variety  of  output  signals.  In  particular,  sinusoidal,  square  and  triangular 
wave  signals  are  generated  by  the  wireless  sensor.  The  signals  measured  by  the  oseilloseope  are 
presented  in  Figure  3.5.  As  ean  be  seen,  the  actuation  interface  was  accurate  (absent  of  noise)  and  could 
be  operated  at  high  sample  rates. 

3.5.  Application  of  IEEE  802.15.4  Protocol  in  an  Agent-Based  Framework 

In  order  to  simplify  the  IEEE  802.15.4  standard,  the  LR-WPAN  architecture  ean  be  broken  down  into  a 
number  of  layers.  Eaeh  layer  corresponds  to  a  specific  part  of  the  standard  and  provides  services  to 
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higher  layers.  The  IEEE  802.15.4  standard  concentrates  only  on  the  bottom  two  layers  of  the 
communication  protocol  stack:  the  physical  (PHY)  and  medium  access  control  (MAC)  layers.  The  PHY 
layer  interacts  directly  with  the  physical  radio  in  order  to  modulate  and  demodulate  data  and  the  MAC 
layer  offers  standardized  packet  structures  for  data  transmission  as  well  as  a  variety  of  methods  for 
autonomous  nodes  to  share  eommon  bandwidth.  The  upper  layers,  consisting  of  a  network  layer  and  an 
application  layer,  are  application  dependent  and  not  defined  explicitly  by  the  standard. 

When  embedded  within  the  operating  system  of  a  wireless  node,  the  upper  layers  represent  any 
application  which  requires  communication  and  data  transfer  services  between  deviees.  In  this  project,  an 
example  application  may  be  control  algorithms  for  naval  systems.  These  upper  layers  are  used  to  aeeess 
the  MAC  layer  for  the  purpose  of  beacon  management,  channel  aeeess,  guaranteed  time  slot  (GTS) 
management,  frame  validation,  and  acknowledged  frame  delivery,  association,  and  disassoeiation.  In  turn, 
the  MAC  layer  utilizes  the  PHY  layer  in  order  to  control  activation  and  deactivation  of  the  radio 
transceiver,  energy  detection,  link  quality  indication,  channel  selection,  and  clear  channel  assessment. 
The  PHY  layer  is  also  used  to  transmit  and  receive  packets  across  the  physical  medium. 

As  modem  naval  systems  become  increasingly  complex,  the  number  of  sensors  and  actuators  in  a  given 
eontrol  system  and  the  subsequent  computational  demands  placed  on  a  centralized  controller  will  grow 
exponentially.  As  a  result,  the  ability  to  control  sueh  a  system  in  a  centralized  manner  is  becoming  ever 
more  difficult.  In  particular,  when  considering  a  wireless  eontrol  system  where  centralized  activity  across 
a  high  density  wireless  network  is  not  at  all  feasible,  an  agent-based  framework  has  been  proposed  in 
which  each  sensing  unit  behaves  in  an  autonomous  manner.  Agents  are  software  abstractions  capable  of 
sensing  their  environment  and  reacting  to  changes  in  that  environment  through  rational  embedded 
intelligence.  Multiple  agent  systems  are  very  attractive  for  controlling  eomplex  systems  because  they 
offer  a  scalable  approach  for  solving  problems  using  decentralized  data. 

In  order  to  aeeommodate  sueh  a  framework  in  a  wireless  setting,  a  layer  of  embedded  firmware  is  needed 
to  serve  as  a  real-time  operating  system  that  allows  upper  software  layers  to  eolleet  and  store  sensor  data, 
issue  actuation  commands,  and  share  state  data  with  other  units  within  the  system.  The  IEEE  802.15.4 
communication  protocol  has  been  effectively  leveraged  for  this  exact  purpose.  In  addition  to  providing  a 
well  defined  application  program  interface  (API)  to  the  upper  layers,  these  embedded  layers  of  the 
communication  protocol  stack  facilitate  communication  between  multiple  agents  aeting  autonomously. 

The  IEEE  802.15.4  communication  protocol  has  many  advantages  that  eater  to  an  autonomous  network. 
To  ensure  that  the  network  is  time  synchronized,  each  unit  is  designed  to  act  in  the  MAC-provided 
beacon-enabled  mode.  Additionally,  channel  access  is  controlled  by  a  slotted  earner  sense  multiple 
access  collision  avoidance  scheme  (CSMA/CA).  In  this  meehanism,  each  node  must  sense  a  free  channel 
at  least  twice  before  transmitting.  The  first  sense  is  delayed  by  a  random  delay,  which  serves  to  reduee 
the  probability  of  collision  when  two  nodes  simultaneously  sense  a  free  ehannel.  When  a  channel  is  busy, 
transmission  docs  not  occur  and  an  additional  random  delay  is  computed.  When  a  packet  does  not  reach 
its  destination,  it  can  be  retransmitted  during  a  new  contention  period. 

The  PHY  and  MAC  layers  have  been  encoded  in  the  C  programming  language  and  loaded  into  the  flash 
memory  of  the  wireless  sensing  unit  microcontroller.  In  the  laboratory,  a  variety  of  MAC  layer  eommand 
and  data  frames  are  transmitted  wirelessly  from  unit  to  unit.  These  command  and  data  frames  are  then 
intercepted  and  interpreted  by  an  IEEE  802.15.4  packet  sniffer  whieh  is  used  to  record  all  paekets  sent 
within  the  network.  These  tests  display  the  functionality  of  the  embedded  MAC  protocol  as  well  as  the 
wireless  connectivity  of  the  system.  A  sereen  shot  of  the  packet  sniffer  results  can  be  seen  in  Figure  3.6. 
In  order  to  further  demonstrate  the  ability  of  the  PHY  layer  to  command  the  physical  radio  to  modulate 
and  demodulate  data,  a  testbed  is  created  whieh  utilizes  this  layer  in  conjunction  with  the  sensing 
capabilities  of  the  wireless  sensing  units.  In  the  laboratory,  a  unidirectional  accelerometer  is  attached  to  a 
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Figure  3.6.  Laboratory  validation  of  full  implementation  of  IEEE  802.15.4  MAC  and  PHY 

layers 


Figure  3.7.  Experimental  setup  and  acceleration  time  history  plot  for  laboratory  validation  of 

PHY  layer  for  sensing  applications. 


vacuum  pump.  The  pump  is  turned  on  and  acceleration  data  is  collected  with  a  wireless  sensing  unit  at  a 
sampling  rate  of  500  Hz.  This  data  is  then  transmitted  back  to  a  centralized  server  and  stored  in  a  local 
data  repository.  Figure  3.7  shows  the  experimental  setup  and  the  time  history  acceleration  record  of  the 
pump  vibration. 
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4.  Model  Updating  for  System  Assessment  and  Damage  Detection: 

In  a  naval  environment,  it  is  imperative  that  distinct  changes  in  the  operational  condition  of  naval  plants 
be  detected,  as  the  detection  of  such  state  changes  will  lead  to  the  initiation  of  market-based  eontrol 
solutions  adopted  to  optimally  allocate  searee  system  resources.  As  such,  this  seetion  outlines  the  ereation 
of  a  novel  data  processing  architecture  used  to  continuously  update  system  models  representing  a  fully- 
funetional  ship  plant  (e.g.,  ehilled  water  supply  system).  Changes  to  sueh  a  plant  are  identified  through 
model  updating  methods  conducted  on  the  wireless  sensor  network,  and  ean  be  indicative  of  damage. 
Model  updating  methods  are  typically  computationally  intense  and  are  often  reserved  for  data  servers 
where  significant  computational  resources  are  located.  In  our  approach,  a  parallel  implementation  of  the 
simulated  annealing  model  update  method  (a  stochastic  search  algorithm)  is  embedded  upon  a  network  of 
Narada  wireless  sensors  to  identify  bifurcations  in  the  system,  perhaps  due  to  battle  damage.  Emphasis  is 
made  to  keep  the  update  time  to  near  real-time  (less  than  ten  minutes)  and  to  ensure  that  the  full 
computational  capabilities  of  the  wireless  sensor  network  (albeit,  highly  distributed)  are  maximized. 

By  viewing  a  wireless  network  as  a  parallel  computer  with  an  unknown  and  possibly  changing  number  of 
processing  nodes,  this  distributed  data  processing  architecture  is  capable  of  performing  complicated  types 
of  data  analysis  while  creating  a  scalable  environment  that  is  not  only  resistant  to  communication  and 
sensor  failure,  but  that  also  becomes  increasingly  efficient  at  higher  nodal  densities.  This  novel 
architecture  functions  by  allowing  a  network  of  sensors  to  autonomously  deteet  and  utilize  the  computing 
resources  of  any  available  wireless  node  on  the  fly.  This  “ad-hoe”  capability  allows  for  increases  in  the 
parallelism  and  efficiency  of  the  architecture  in  real-time,  and  ean  be  used  to  reform  or  “self-heal”  the 
network  in  the  wake  of  any  communication  and/or  sensor  failures. 

In  order  to  examine  the  data  processing  capabilities  of  this  novel  architecture  from  an  damage  detection 
perspective,  a  parallelized  version  of  the  simulated  annealing  (SA)  stoehastic  optimization  method  is 
designed  for  implementation  within  a  distributed  WSN.  One  of  the  reasons  that  the  SA  algorithm  is 
ehosen  for  parallelization  is  that  it  can  be  applied  to  many  of  the  optimization  problems  that  may  arise  in 
the  naval  environment.  In  this  study,  the  wireless  parallel  SA  (WPSA)  method  is  used  within  a  WSN  for 
the  updating  of  an  analytical  model  of  a  shipboard  ehilled  water  system. 

4.1.  Background  on  combinatorial  optimization  by  simulated  annealing: 

One  of  the  most  studied  areas  in  computational  engineering  is  that  of  combinatorial  optimization  (CO). 
This  field  involves  developing  efficient  methods  for  finding  the  maximum  or  minimum  value  of  any 
function  with  a  large  number  of  independent  variables.  CO  problems  are  typieally  very  difficult  to  solve 
computationally,  as  an  exact  solution  often  requires  a  number  of  computational  steps  that  grows  faster 
than  any  finite  power  of  the  size  of  the  problem.  As  such,  it  is  often  desirable  in  engineering  applications 
to  quiekly  find  good  approximations  to  the  optimal  solution  instead  of  expending  the  time  and  resources 
required  to  find  an  absolute  global  optimum.  Unfortunately,  even  approximate  solutions  ean  sometimes 
be  difficult  to  find,  as  most  relevant  search  strategies  involve  iterative  improvement,  and  as  sueh,  have  a 
tendency  to  get  stuck  in  loeal  (not  global)  optima.  However,  in  the  1980’s,  several  algorithms  derived 
from  physical  and  biological  systems  were  developed  for  finding  near-global  optima  in  functions 
containing  many  local  optima  (Bounds  1987).  One  of  these  methods  is  the  simulated  annealing  (SA) 
optimization  technique,  first  presented  by  Kirkpatrick  et  ai  (1983). 

SA  was  developed  out  of  the  observation  that  a  connection  eould  be  made  between  CO  and  the  behavior 
of  physical  material  systems  in  thermal  equilibrium  at  a  finite  temperature.  In  material  physics, 
experiments  that  determine  the  low-temperature  state  of  a  material  are  performed  by  first  melting  the 
substanee,  and  then  slowly  lowering  the  substanee’s  temperature,  eventually  spending  a  long  time  at 
temperatures  near  freezing.  This  annealing  procedure  allows  the  substance  to  eventually  obtain  an  optimal 
thermal  energy  state  amongst  an  almost  infinite  number  of  possible  atomistic  configurations.  Assuming 
that  a  method  exists  for  determining  the  energy  of  a  physical  system  in  a  specific  atomistic  configuration, 
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this  physical  annealing  procedure  can  be  viewed  as  a  CO  problem  where  the  objective  is  to  find  the 
globally  minimal  energy  state  of  the  material’s  atoms. 


As  such,  by  borrowing  ideas  from  the  natural  annealing  process,  a  “simulated’’  version  of  the  annealing 
method  can  be  developed  to  quickly  obtain  good  approximate  solutions  to  CO  problems  where  the 
objective  is  to  find  a  globally  minimal  value  of  some  optimization  function.  This  is  done  by  viewing  the 
value  of  the  function  to  be  optimized  as  the  physical  system  “energy”,  introducing  an  “effective” 
annealing  temperature  which  will  simulate  the  material  cooling  process,  and  utilizing  the  Metropolis 
procedure  (explained  below')  to  avoid  premature  convergence  on  local  optima,  which  is  the  key  to  the 
effectiveness  of  the  generalized  annealing  process. 

In  1953,  Metropolis  et  al.  created  an  algorithm  that  can  probabilistically  simulate  a  collection  of  atoms 
converging  on  thermal  equilibrium  at  a  set  temperature.  At  each  step  in  this  algorithm,  a  randomly 
selected  atom  is  displaced  a  small,  random  distance,  and  the  resulting  change  in  system  energy  (AE)  is 
computed.  If  A E  <  0,  this  disturbance  is  accepted.  Otherwise,  if  A E  >  0,  the  new'  configuration  will  be 
accepted  with  the  following  probability: 


Pr (accept)  ~  e 


-A £ 

k.T 


(4.1) 


where  T  is  the  temperature  of  the  system  and  is  Boltzmann’s  constant.  If  the  new  configuration  is 
accepted,  the  next  step  of  the  search  continues  with  that  atom  displaced.  Otherwise,  if  the  new 
configuration  is  not  accepted,  the  next  step  in  the  search  continues  using  the  original  atomistic 
configuration.  By  repeating  this  procedure  many  times.  Metropolis  simulates  the  thermal  motion  of  atoms 
subjected  to  a  constant  temperature,  and  mimics  the  probabilistic  process  by  which  nature  avoids 
premature  convergence  on  suboptimal  configurations. 

As  proposed  by  Kirkpatrick  et  al.  (1983),  “simulated”  annealing  can  be  used  in  the  context  of  CO  by 
representing  each  possible  configuration  of  optimization  function  parameters  as  a  distinct  state,  s.  The 
objective  of  the  annealing  process  is  to  find  a  system  state  that  minimizes  the  value  of  an  optimization 
function,  E(s).  In  order  to  help  avoid  convergence  on  a  sub-optimal  minimum,  the  Metropolis  framework 
can  be  applied  to  the  SA  procedure  by  generating  a  new  state,  snew,  by  altering  the  value  of  one  function 
parameter  at  random.  The  objective  function  value  of  this  new'  state,  E(sneJ),  is  then  compared  with  the 
objective  function  value  of  the  old  state,  E(saiJ),  and  the  newr  state  is  probabilistically  accepted  or  rejected 
based  on  the  criterion  presented  in  Equation  4.1.  When  SA  is  implemented  within  a  computing  machine, 
the  probability  of  a  new  system  state  being  accepted  at  a  given  temperature  can  be  stated  as  follows: 
accept  a  new  state,  snewy  if  and  only  if: 


)<E{snld)+T.\\n(u] 


(4.2) 


where  U  is  a  uniformly  distributed  random  variable  between  0  and  1.  The  addition  of  the  7’*|ln((7)|  term 
allows  the  system  to  periodically  accept  a  sub-optimal  state  in  hopes  of  avoiding  premature  convergence 
on  a  local  optima.  A  standard  SA  cooling  schedule  begins  the  optimization  process  by  assigning  a  high 
initial  temperature  To  and  then  letting  the  Metropolis  algorithm  run  for  A  iterations.  During  each  iteration, 
a  new  psuedorandomly  generated  state  is  created  by  modifying  one  of  the  optimization  parameters,  and 
the  newly  generated  state  is  either  rejected  or  accepted  based  on  the  Metropolis  criterion  (Equation  4.2). 
After  A  iterations,  the  temperature  of  the  system  is  reduced  by  a  factor  of  p ,  such  that  Tnew  -  p'To}(j,  and  A 
additional  iterations  will  be  run  at  the  new,  lower  temperature  (Tncii).  This  process  continues  until  the 
temperature  is  sufficiently  low  that  very  few'  new  states  are  accepted,  meaning  that  a  globally  optimal 
state  has  likely  been  found  and  the  system  has,  in  essenee,  frozen. 
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Figure  4.1:  Random  state  generation  for  a  two-dimensional  search  problem  using 
(a)  standard  SA,  (b)  BSA,  and  (c)  WPSA  (Zimmerman  and  Lynch  in  press). 


Since  Kirkpatrick  et  ah  first  published  the  SA  methodology  in  1983,  countless  variations  on  the  original 
algorithm  have  been  seen  in  the  literature.  For  each  specific  optimization  problem,  it  seems,  a  different 
variant  on  the  traditional  SA  method  provides  the  quickest  convergence  and  the  most  accurate  results.  As 
such,  it  is  important  to  note  that  the  WPSA  methodology  proposed  herein  for  use  in  wireless  sensor 
networks  can  be  effectively  utilized  in  conjunction  with  almost  any  variant  on  the  SA  method.  However, 
for  the  model  updating  problem  studied  in  this  paper,  a  modification  on  the  blended  simulated  annealing 
(BSA)  algorithm  proposed  by  Levin  and  Lieven  (1998)  is  exclusively  utilized.  The  BSA  algorithm 
deviates  from  the  standard  SA  methodology  in  the  way  in  which  it  creates  randomly  generated  states.  In 
standard  SA,  new  states  are  generated  by  randomly  choosing  one  annealing  parameter  and  assigning  it  a 
new  value  chosen  uniformly  from  within  the  parameter’s  valid  range  (Figure  4.1(a)).  In  the  BSA 
algorithm,  however,  this  standard  type  of  state  generation  is  alternated  every  other  step  with  a  ‘’radius 
adjustment”  approach,  where  all  annealing  parameters  are  changed  by  choosing  a  random  point  on  a 
hypersphere  that  is  a  fixed  radius  away  from  the  previous  annealing  state  (Figure  4.1(b)).  This  method 
requires  two  separate  annealing  temperatures,  one  for  the  standard  SA  adjustment  and  one  for  the  radius 
adjustment.  For  this  study,  the  BSA  algorithm  is  modified  slightly  such  that  instead  of  choosing  a  point 
that  lies  on  a  fixed  radius  from  a  previous  annealing  state,  all  annealing  parameters  are  randomly  assigned 
new  values  that  reside  within  a  given  radius  from  the  individual  parameter’s  current  assignment  (Figure 
4.1(c)).  Then,  the  radius  itself  is  treated  as  a  variable  in  the  SA  process  much  like  the  annealing 
temperature.  It  starts  w  ith  a  high  value  near  1.0  (such  that  the  entirety  of  each  parameter’s  valid  range  can 
be  searched),  and  as  time  progresses,  the  searchable  radius  is  reduced  such  that  the  SA  search  focuses 
increasingly  on  values  that  are  close  to  the  currently  optimal  state.  This  improves  upon  the  BSA  algorithm 
by  eliminating  the  w  asteful  interrogation  of  search  states  far  away  from  the  currently  optimal,  especially 
later  in  the  search  as  a  final,  optimal  solution  is  converged  upon. 
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4.2.  Wireless  parallel  simulated  annealing: 

When  considering  performing  CO  tasks  on  a  wireless  sensor  network,  SA  may  at  first  appear  to  be  an 
excellent  candidate  for  a  stochastic  search  procedure.  Because  a  search  using  SA  requires  only  a 
negligible  two  or  possibly  three  states  to  be  stored  in  memory  at  any  one  time,  SA  is  extremely  attractive 
in  the  wireless  setting  where  memory  capacity  within  most  prototypes  is  limited.  However,  the 
computational  costs  of  implementing  SA,  which  may  require  a  value  of  E  to  be  determined  at  hundreds  of 
thousands  of  randomly  generated  states  in  order  to  converge  on  an  optimal  solution,  can  be  staggering. 
When  implemented  within  a  single  wireless  sensing  device,  where  processing  speed  is  usually  only  a 
fraction  of  that  of  an  ordinary  personal  computer,  this  is  a  potentially  debilitating  problem. 

In  order  to  mitigate  the  computational  demands  imposed  by  SA,  many  researchers  have  developed 
parallel  SA  techniques  that,  when  run  on  a  large  number  of  processors,  can  successfully  increase  the 
speed  with  which  a  solution  to  a  CO  problem  can  be  obtained  (Greening  1990).  However,  most  of  these 
methods  require  communication  between  processors  both  before  and  after  each  random  state  is  generated. 
In  the  wireless  setting,  where  battery  preservation  is  a  high  priority  and  communication  bandwidth  is 
limited,  this  type  of  constant  communication  negates  the  advantages  of  parallelism  and  represents  a  poor 
use  of  battery  power.  In  this  study,  a  parallel  SA  procedure  is  created  that  utilizes  the  computational 
resources  distributed  across  large  wireless  sensing  networks  while  minimizing  the  communication 
demands  of  the  parallel  algorithm.  This  is  done  by  taking  advantage  of  the  fact  that  the  SA  process 
typically  rejects  more  states  than  it  accepts,  especially  as  the  annealing  temperature  is  lowered  and  the 
algorithm  converges  on  a  solution.  Specifically,  the  traditionally  serial  SA  search  problem  (which  is 
continuous  across  all  temperature  steps)  can  be  broken  into  a  set  of  smaller  search  trees,  each  of  which 
corresponds  to  a  given  temperature  step  and  begins  with  the  globally  optimal  state  assignment  so  far 
detected  at  the  preceding  temperature  step.  Each  smaller  search  problem  can  then  be  assigned  individually 
to  any  available  sensor  in  the  network,  and  thus  multiple  temperature  steps  can  be  searched  concurrently. 
This  concept  is  displayed  graphically  in  Figure  4.2.  One  of  the  great  advantages  of  this  methodology  is 
that,  given  the  ad-hoc  communication  capabilities  of  many  wireless  sensing  devices,  these  individualized 
search  trees  can  be  distributed  in  real-time  to  any  available  processor  within  the  sensing  network.  Because 
the  ad-hoc  assignment  and  reassignment  of  search  problems  can  allow  for  individual  nodes  to  drop  from. 


Figure  4.2:  A  simple  serial  SA  seareh  tree,  shown  up  to  the  fourth  temperature  step  and  its 
corresponding  YVPSA  seareh  trees,  assigned  to  w  ireless  sensors  (Zimmerman  and  Lynch  in  press). 
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or  appear  in  the  network  mid-scarch,  this  parallelized  updating  method  is  incredibly  valuable  in  systems 
where  sensor  or  communication  reliability  may  be  in  question. 

4.2.1.  Wireless  implementation  of  the  WPSA  algorithm: 

In  the  parallel  SA  implementation  used  in  this  study  (WPSA),  a  computational  task  requiring  SA 
optimization  is  first  assigned  to  any  one  available  sensing  unit,  along  with  a  user-defined  initial 
temperature,  Tq.  This  first  wireless  sensor,  n0,  then  beacons  the  network,  searching  for  other  sensors 
available  for  data  processing.  If  a  seeond  sensing  node,  nu  is  found,  the  first  sensor,  n(h  w  ill  assign  the  S  A 
search  tree  starting  at  the  next  temperature  step,  Tu  to  the  second  sensor,  n y,  passing  along  its  eurrent 
information  regarding  the  most  optimal  system  state  yet  visited.  This  process  continues  until  no  sensors 
remain  available  for  data  processing. 

If  a  given  sensor,  deteets  an  optimal  solution,  (i.e.,  no  newr  states  are  accepted  at  the  temperature  step 
sensor  nt  is  investigating),  it  will  order  the  rest  of  the  network  to  discontinue  the  SA  search,  and  w  ill  alert 
the  network  end-user  of  the  discovered  results.  However,  if  sensor  n{  finishes  its  part  of  the  SA  seareh 
without  having  converged  on  a  solution  (i.e.,  new  states  are  still  being  accepted),  it  will  alert  its  sueeessor, 
/7/+/,  that  no  solution  was  found  at  temperature  step  7),  and  sensor  m  will  again  make  itself  available  to  the 
network  for  computation  on  a  lower  temperature  step.  While  WPSA  functions  autonomously  without 
need  for  a  centralized  controller,  the  WSU  assigned  to  the  highest  temperature  step  at  any  given  time 
keeps  track  of  seareh  progress  and  alerts  the  user  when  the  search  has  been  completed.  Beeause  of  the 
self-healing  capabilities  of  many  WSNs,  this  parallel  algorithm  will  always  adapt  in  order  to  utilize  the 
maximum  number  of  processing  nodes  available  at  any  one  time,  even  if  some  sensors  drop  in  and  out  of 
the  network  during  computation. 

As  the  WPSA  search  continues,  information  regarding  newly  found,  increasingly  optimal  states  is 
disseminated  downwards  through  the  network,  sueh  that  all  sensors  are  eognizant  of  any  seareh  progress 
that  has  been  made  at  higher  temperature  steps.  This  allows  all  sensors  to  maximize  the  effectiveness  of 
their  search  at  a  given  temperature  step,  and  maintains  the  continuity  of  the  serial  SA  process.  Spceifieally, 
when  a  sensor  deteets  a  state,  s ,  with  a  lower  optimization  function  value  than  that  of  any  other  known 
state,  it  will  immediately  propagate  this  information  downward  to  the  sensor  direetly  below  it  in  the 
seareh  tree  (its  ehild).  If  the  propagated  state  information  also  represents  the  minimal  value  of  the 
objective  function  that  the  ehild  has  found  so  far,  the  ehild  will  then  restart  its  N  seareh  iterations  from  the 
newly  found  minimum  state  and  inform  the  sensor  direetly  belowr  it  of  this  newly  discovered  state. 
However,  if  a  ehild  receives  a  state,  sp,  from  a  parent,  and  the  ehild  has  already  randomly  generated  a 
state,  sc ,  that  yields  a  lower  objective  function  value  than  sf„  (E(sc)  <  E(sp))r  that  child  will  merely  restart 
its  SA  iterations  given  its  current  search  state,  sc ,  without  passing  any  information  on  to  its  successor.  In 
this  way,  it  is  assured  that  each  temperature  step  is  thoroughly  searched  given  the  complete  information 
obtained  at  the  preceding  temperature  step.  While  this  does  result  in  an  increase  in  the  total  number  of  SA 
iterations  required  to  reach  a  solution  over  the  serial  SA  procedure,  the  additional  randomly  generated 
states  at  many  (if  not  all)  temperature  steps  slightly  increases  the  probability  that  a  “better”  solution  will 
be  found  than  otherw  ise  possible. 

4.2.2.  Illustrative  example  of  the  WPSA  algorithm: 

Figure  4.3  illustrates  the  distribution  of  one  example  parallelized  simulated  annealing  task  over  a  network 
of  four  wireless  sensing  units.  This  task  has  an  initial  global  minimum  objective  function  value  of  E  =  10, 
and  is  assigned  by  the  user  to  wireless  sensing  unit  1  (WSU  1).  Simultaneous  to  this  assignment,  the  user 
also  alerts  all  other  sensors  that  they  should  make  themselves  available  for  computation.  After  reecipt  of 
this  task  assignment,  WSU  1  recognizes  that  WSU  2  is  an  available  computational  node,  and  orders  this 
unit  to  perform  N  seareh  iterations  at  the  second  temperature  step,  starting  with  WSU  1  ’s  current  global 
state  (with  a  global  minimum  of  E  =  10).  In  a  similar  way,  WSU  2  assigns  the  third  temperature  step  to 
WSU  3,  and  WSU  3  assigns  the  fourth  temperature  step  to  WSU  4. 
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After  searching  approximately  Nl 4  SA-generated  search  states,  WSU  3  detects  a  state  with  an  objective 
function  value  of  8.  It  immediately  passes  this  state  information  along  to  its  child,  WSU  4.  Because  WSU 
4  has  a  current  global  minimum  value  of  E=  10,  WSU  4  restarts  its  search  of  N  SA-generated  states  at  the 
fourth  temperature  step  with  this  updated  information.  WSU  4  has  no  children,  so  the  propagation  of  this 
new  state  stops  when  it  reaches  WSU  4. 

Soon  thereafter,  WSU  1  detects  a  SA-generated  state  with  an  objective  function  value  of  E  =  9.  This  is 
lower  than  its  current  minimum  value  of  E  =  10,  so  it  informs  its  child,  WSU  2  of  the  newly  found  state. 
WSU  2  recognizes  E  -  9  as  a  new  global  minimum,  so  it  restarts  its  search  at  the  second  temperature  step 
with  this  information  and  passes  the  updated  values  along  to  its  child  (WSU  3).  WSU  3,  however,  has 
already  detected  a  global  minimum  of  E  =  8,  and  thus  it  simply  restarts  its  search  of  N  SA-generated 
states  with  its  current  state  information  (without  needing  to  send  any  updated  information  to  its  child). 

When  WSU  1  finishes  its  search  of  N  SA-generated  states,  it  alerts  its  child  (WSU  2)  that  it  has  not  found 
a  globally  optimal  solution,  and  it  disengages  from  the  search  process.  At  this  point,  however,  WSU  1 
broadcasts  its  availability  to  the  other  nodes  in  the  network.  WSU  4,  which  is  in  need  of  a  child  node, 
assigns  the  SA  task  at  the  fifth  temperature  step  to  WSU  1,  given  its  current  state  (with  a  global  minimum 
value  of  E-  8). 

This  process  continues  for  several  more  temperature  steps  until  WSU  1  detects  a  globally  optimal  state 
(i.e.  it  finds  a  state  with  an  objective  function  value  of  E-  0).  At  this  point,  WSU  1  broadcasts  its  find  to 
the  network,  thereby  stopping  all  other  computation,  and  it  alerts  the  user  that  a  globally  optimal  state  has 
been  found. 

4,3,  Experimental  Validation  of  the  VVPSA  algorithm  within  a  Naval  Chilled  Water  System: 

To  validate  the  performance  of  the  proposed  distributed  WPS  A  process  for  damage  detection,  a  set  of 
experiments  were  run  on  simulated  data  from  a  very  simple  Navy  shipboard  chilled  water  system  (shown 
in  Figure  4.4).  The  pipe  network  in  this  system  contains  sixteen  pipes  of  equal  diameter  but  varying 
length.  The  system  is  driven  by  two  identical  pumps,  located  along  pipes  three  and  four,  and  there  are 
four  flowmeters  measuring  flow  within  the  network.  Each  of  these  flowmeters  is  monitored  by  a  single 
wireless  sensing  unit. 

An  analytical  model  of  the  pipe  network  described  above  is  designed  and  embedded  within  each  sensing 
node.  Given  a  particular  system  configuration,  this  software  can  determine  the  expected  flowrate  and  head 
loss  in  each  pipe.  By  conserving  mass  at  each  pipe  junction  and  using  the  linearized  Darey-Weisback 
equation  to  relate  head  loss  in  loops,  a  set  of  linear  equations  describing  the  piping  network  can  be 
automatically  formulated.  These  equations  can  then  be  solved  using  Gauss-Jordan  elimination,  and  values 
for  flowrate  can  be  found.  However,  because  head  losses  in  this  environment  are  dependent  on  flowrates 
and  flowrates  are  dependent  on  head  losses,  this  process  must  iterate  until  a  solution  is  converged  upon. 
Once  embedded  within  the  Narada’s  Atmegal28  processor,  each  execution  of  the  Gauss-Jordan  process  is 
timed  and  found  to  take  roughly  0.73  seconds. 
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Annealing  Temperature  and  Annealing  Temperature  and 

Network  Topology  Task  Distribution 


Figure  4.3:  One  parallel  simulated  annealing  task  running  on  four  wireless  sensors 
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Figure  4.4:  A  simple  Navy  shipboard  chilled  water  system. 

Because  flowrate  data  is  constantly  available  at  several  locations  throughout  the  pipe  network,  an 
objective  function  can  easily  be  created  to  compare  experimental  observations  with  analytical  model 
projections: 


£=IV(e--&'„J 

i 


(4.3) 


where  E  is  the  objective  difference  between  analytical  and  experimental  results  and  Q  anal  and  Q sense  arc 
the  analytical  and  sensed  flowrates,  respectively,  in  pipe  /. 

Using  the  SA  method,  this  objective  function  can  be  repeatedly  evaluated  given  different  values  for  a  set 
of  model  properties  until  an  “updated”  configuration  of  parameters  is  found  that  minimizes  the  objective 
difference  between  expected  and  observed  flowrates.  Because  changes  in  this  objective  function  can  be 
interpreted  as  a  direct  indicator  of  damage,  the  choice  of  updating  parameters  is  extremely  important  in 
terms  of  the  ability  to  detect  and  localize  system  bifurcations.  In  the  case  of  a  chilled  water  system,  it  is 
advantageous  to  choose  updating  parameters  that  strongly  influence  flowrate,  such  as  pipe  diameter.  By 
doing  so,  it  becomes  possible  to  look  for  changes  in  pipe  properties  as  an  indication  of  suboptimal 
performance. 

For  the  sake  of  simplicity,  it  is  assumed  that  all  damage  in  this  system  is  local  to  sector  three,  which 
contains  pipes  12,13,14,  and  15.  In  order  to  account  for  a  reduction  in  flow  rate  that  may  result  from  a 
blocked  or  leaking  pipe,  pipe  diameters  for  each  of  these  four  pipes  are  designated  as  the  updating 
parameters  for  this  problem.  All  pipes  in  the  system  arc  initialized  with  a  diameter  of  0.5  inches.  In  order 
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to  introduce  simulated  damage  to  the  system,  sensor  results  are  simulated  with  pipe  13  having  an  80% 
reduction  in  pipe  diameter. 

Having  obtained  sensor  readings  for  the  damaged  state,  a  central  server  sets  up  the  model  updating 
problem  and  chooses  one  sensing  unit  to  serve  as  the  initial  network  coordinator.  This  unit  then  searches 
for  an  available  node  to  assign  additional  tasks  to,  and  the  SA  process  continues  in  an  ad-hoc  fashion  until 
a  set  of  updated  parameters  arc  converged  upon.  A  time  series  representation  of  pipe  diameter 
convergence  can  be  seen  in  Figure  4.5. 

In  order  to  evaluate  the  effectiveness  of  the  parallelized  updating  algorithm,  an  entire  model  update  is 
completed  several  times  with  varying  numbers  of  Narada  nodes  made  available  for  computation.  It  can  be 
seen  from  Figure  4.6  that  this  type  of  parallelization  can  drastically  decrease  the  time  required  to  find 
updated  model  properties  and  detect  damage  within  a  system.  It  can  also  be  seen  from  this  figure  that 
increasing  the  number  of  computing  sites  also  improves  the  accuracy  of  the  final  solution.  Figure  4.7 
shows  the  ability  of  this  algorithm  to  accurately  detect  and  locate  system  bifurcations.  It  can  be  seen  that 
the  SA  method  rightly  identified  a  significant  diameter  reduction  in  pipe  13,  while  pipes  12,14,  and  15 
only  observe  minor  variations. 
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Figure  4.5:  SA  pipe  diameter  convergence. 
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Relative  Time  to  Reach  Converged  Solution 


Number  of  Wireless  Sensing  Units  Number  of  Wireless  Sensing  Units 

Figure  4.6:  Speed  and  accuracy  improvements  using  WPSA. 
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Figure  4.7:  Speed  and  accuracy  improvements  using  WPSA. 
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5.  Experimental  Testing  of  Narada  Wireless  Sensors  in  a  Chilled  Water  Plant: 

In  order  to  experimentally  verify  the  ability  of  the  Narada  wireless  sensor  to  perform  sensing,  aetuation 
and  control  tasks  relevant  to  the  Naval  environment,  the  ONR  tabletop  chilled  water  demonstration 
system  was  chosen  as  an  appropriate  testbed.  Field  validation  of  the  Narada  sensing  platform  was 
performed  in  two  phases,  with  Phase  I  focusing  on  the  integration  of  the  wireless  sensing  environment 
within  the  chilled  water  system  and  Phase  11  investigating  the  ability  of  a  network  of  Narada  wireless 
sensors  to  detect  and  locate  plant  damage  using  the  model  updating  methodology  presented  in  Section  4. 

5.1.  Phase  I  (Proof  of  Concept)  Field  Work  Statement: 

This  section  of  the  report  details  the  results  of  field  testing  performed  on  June  25th  and  26th,  2008,  at  the 
Naval  Surface  Warfare  Center  in  Philadelphia,  PA.  On  these  dates,  the  basic  interface  capabilities  of  the 
University  of  Michigan’s  Narada  wireless  sensing  system  were  tested  using  the  ONR  tabletop  chilled 
water  demonstration  system.  The  purpose  of  these  field  tests  was  to  validate  the  ability  of  the  Narada 
wireless  system  to  record  flow  data  within  a  chilled  water  system,  to  archive  this  data  for  later  access,  and 
to  lay  the  groundwork  for  the  use  of  embedded  model  updating  techniques  to  detect  and  locate  pipe  or 
pump  degradation  within  the  chilled  water  system. 

5.1.1.  Phase  1  Objectives: 

The  objectives  of  the  Phase  I  field  work  can  be  defined  as  follows: 

1.  Interface  Narada  wireless  sensors  with  the  flowmeters  already  installed  on  the  ONR  tabletop 
chilled  water  demonstration  system. 

2.  Validate  the  ability  of  the  Narada  sensors  to  collect  analog  flowmeter  data  from  the  tabletop 
demonstrator  and  convert  it  into  digital  flow  readings. 

3.  Collect  full-system  flow  data  from  the  tabletop  demonstrator  in  an  undamaged  state. 

4.  Introduce  simulated  damage  to  the  tabletop  demonstrator  through  the  use  of  rupture  valves,  pump 
shorts,  and  valve  closures. 

5.  Collect  full-system  flow  data  from  the  tabletop  demonstrator  in  a  damaged  state. 

5.1.2.  Integration  of  Narada  wireless  sensors  with  tabletop  flowmeters: 

The  ONR  tabletop  chilled  water  demonstration  system  (shown  in  Figure  5.1)  was  developed  as  a  testbed 
on  which  to  validate  new  sensing  and  control  technologies.  This  tabletop  system,  which  emulates  many  of 
the  electrical  and  mechanical  aspects  of  a  shipboard  chilled  water  system,  contains  four  pumps  which 
provide  chilled  water  through  a  series  of  42  pipes  to  two  vital  and  four  non-vital  “loads”.  In  reality,  these 
loads  are  things  vital  to  shipboard  functionality  (such  as  electrical  or  mechanical  equipment),  but  on  the 
demonstrator  they  are  simply  visual  flow  indicators.  Each  pipe  that  provides  flow  to  a  pump  or  a  load  is 
instrumented  with  a  flowmeter  (for  a  total  of  10  sensors),  and  a  set  of  26  valves  arc  in  place  to  control  the 
flow  to  different  sectors  of  the  system.  Two  rupture  valves  are  also  installed  in  order  to  create  a 
controllable  way  of  simulating  damage  within  the  pipe  network.  The  testbed  also  contains  an  electrical 
system  built  to  emulate  shipboard  conditions,  but  for  the  purposes  of  this  study  only  the  mechanical  pipe 
network  is  of  interest.  Details  of  the  equipment  used  to  interface  with  the  Narada  sensing  system  can  be 
seen  below: 

•  DigiFlow  DFS-2W  Flowmeter  (Figure  5.2) 

o  Excitation  Voltage:  5V 

o  Measurement  Range  with  Narada  Sensor:  0  mL/sec  -  27.7  mL/sec 

•  Greylor  PQ-12  DC  Gear  Pump  (Figure  5.3) 

o  Characteristic  Equation:  hp  =  -268.527Q2  +  86.78Q  +  40.272 
■  Where  0  is  in  GPM  and  hp  is  in  ft 

•  STC  2W040-3/8-D  Solenoid  Valve  (Figure  5.4) 

o  Excitation  Voltage:  12-24V 


Page  25  of  61 


Figure  5.1.  ONR  tabletop  ehilled  water  demonstration  system,  (a)  plan  view  and  (b)  schematic. 


Figure  5.2.  Digiflow  DFS-2W  Flowmeter,  modified  for  connection  to  Narada  wireless  sensor. 
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Figure  5.3.  Greylor  PQ-12  DC  Gear  Pump. 


Figure  5.4.  STC  2W040-3/8-D  Solenoid  Valve. 


Each  of  the  flowmeters  in  the  tabletop  system  functions  in  the  same  way:  as  fluid  flows  through  a  meter, 
it  first  passes  a  fixed  worm.  This  worm  imparts  a  spiral  flow  which,  in  turn,  spins  a  rotor  on  a  virtually 
friction-free  bearing.  As  they  spin,  the  rotor  blades  interrupt  an  infrared  beam  eoming  from  a  fixed 
external  electronics  mounting  assembly,  thus  generating  a  square  wave  output  signal.  Each  flowmeter  is 
powered  by  a  1 2V  DC  power  supply.  However,  through  the  use  of  an  external  resistor,  the  nominal  output 
voltage  of  the  sensor  is  approximately  3.5V.  Each  time  the  rotor  blades  spin,  an  open  circuit  is  ereated 
and  this  sensor  voltage  drops  briefly  to  OV.  An  example  of  this  analog  output  data  can  be  seen  in  Figure 
5.5.  One  Narada  wireless  sensor  is  utilized  to  monitor  each  flowmeter  (Figure  5.6).  During  the  course  of 
this  field  test,  it  was  found  that  the  voltage  supply  being  utilized  was  providing  an  excessive  amount  of 
noise  to  the  sensor  readings.  This  noise  was  likely  present  because  the  original  power  supply  was  being 
utilized  to  power  not  only  the  flow  meters,  but  all  of  the  demonstrator’s  valves  and  pumps.  The  resulting 
noise  prevented  the  square  wave  signal  from  reaching  the  OV  level,  and  made  it  very  difficult  to  extract 
digital  flowrate  information  from  the  analog  sensor.  An  example  of  this  noisy  output  can  be  found  in 
Figure  5.7.  A  different  power  supply  was  brought  in  for  the  sole  purpose  of  powering  the  system's  ten 
flowmeters,  and  this  solution  appeared  to  solve  the  problem  and  produce  clean  square  wave  output 
(Figure  5.5). 
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In  order  to  convert  this  pulsing  analog  voltage  data  into  flowrate,  specialized  data  acquisition  code  is 
written  for  the  Narada  wireless  sensors.  When  asked  to  record  flow  data,  the  sensors  directly  sample  the 
flowmeter  signal  at  a  very  high  sampling  rate  (5000Hz),  counting  how  frequently  a  pulse  occurs  by 
looking  for  “falling  edges”  or  pairs  of  data  points  where  the  first  is  above  3  V  and  the  second  is  below  3V. 
The  amount  of  time  between  pulses  is  determined  and  a  conversion  factor  is  used  to  store  a  “current 
flowrate”  value  in  memory  (Equation  5.1). 
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Then,  this  “current  flowrate”  value  can  be  sampled  much  more  slowly  by  the  traditional  data  acquisition 
algorithm.  This  method  allows  the  user  to  maximize  the  sampling  speed  of  the  Narada  hardware  while  not 
needlessly  wasting  communication  bandwidth  and  hard  drive  space  by  storing  large  amount  of  extremely 
high  frequency  data.  Lastly,  in  order  to  account  for  variability  and  resolution  error  in  the  edge  counting 
algorithm,  the  converted  flowrate  is  smoothed  using  a  rolling  average.  A  visual  example  of  the  Narada \s 
final  flowTate  output  (sampled  at  200  Hz)  can  be  seen  in  Figure  5.8.  The  blue  dots  in  this  figure  represent 
the  converted  and  sampled  flowrates  and  the  red  line  represents  the  rolling  average. 


Flowmeter  Voltage  Data  -  WSU  9 


F  igure  5.5.  Sample  flowmeter  voltage  data,  collected  with  a  Narada  wireless  sensor. 
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Figure  5.6.  Narada  wireless  sensors  set  to  monitor  tabletop  flowrates. 


Flowmeter  Voltage  Data  -  WSU56 


Figure  5.7.  Sample  noisy  flowmeter  voltage  data,  collected  with  a  Narada  wireless  sensor. 
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Converted  Flowmeter  Data  -  WSU28 
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Figure  5.8.  Sample  converted  flowmeter  data,  collected  and  converted  with  a  Narada  wireless 

sensor. 


Because  the  rotor  blades  spin  very  fast  (60,000  pulses  per  liter  of  passing  water),  a  very  high-frequency 
square-wave  output  signal  is  created,  even  at  the  relatively  lowr  flowrates  seen  on  the  demonstrator 
(typical  5-25  mL/see  flowrates  create  square-wave  frequencies  from  300- 1500Hz).  When  these  high- 
frequcney  signals  arc  compounded  with  the  fact  that  the  square-wave  has  a  very  low'  duty  cycle  (<20%),  it 
beeomes  very  difficult  to  capture  the  voltage  drops  that  oeeur  as  the  rotor  blades  spin.  In  fact,  while  the 
system  is  capable  of  measuring  flowrates  up  to  27.7  mL/see,  it  was  found  that  flow'  rates  above  22  mL/see 
were  often  improperly  recorded  when  using  the  Narada  system,  even  when  sampled  at  the  wireless 
sensor's  maximum  speed  (10,000  Hz).  An  example  Narada  output  from  a  high  flowrate  pipe  is  shown  in 
Figure  5.9.  It  ean  be  seen  that  very  few'  voltage  drops  are  recorded,  when  contrasted  with  the  lower 
flowrate  output  in  Figure  5.5.  It  is  unknown  whether  the  underlying  problem  lies  in  insufficient  sampling 
speeds,  or  whether  the  flowmeters  themselves  are  simply  not  capable  of  reading  flowrates  up  to  their 
documented  maximum  (75  mL/see).  If  it  is  a  sampling  issue,  this  problem  eould  be  remedied  in  the 
future  on  the  Narada  system  by  utilizing  a  faster  timing  crystal  and/or  a  faster  microcontroller  to 
command  the  high-speed  analog  to  digital  conversion,  by  utilizing  an  additional  microcontroller  for  the 
sole  purpose  of  converting  from  analog  to  flow  data,  or  by  selecting  a  different  flowmeter. 
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Excessive  Flowrate  Flowmeter  Data  -  WSU  29 
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Figure  5.9.  Excessive  flowrate  data  (>25mL/sec)  sampled  at  10,000  Hz. 


5.2  Phase  II  (Model  Updating  and  Plant  Damage)  Field  Work  Statement: 

This  section  of  the  report  details  the  results  of  field  testing  performed  on  December  18th  and  19th,  2008,  at 
the  Naval  Surface  Warfare  Center  in  Philadelphia,  PA.  On  these  dates,  the  model  updating  and  damage 
detection  capabilities  of  the  University  of  Michigan’s  Narada  wireless  sensing  system  were  tested  using 
the  ONR  tabletop  chilled  water  demonstration  system.  The  purpose  of  these  field  tests  was  to  record  a 
large  amount  of  flow  data  from  within  a  chilled  water  system  when  that  system  is  subjected  to  several 
types  of  rupture  and  when  that  system  is  being  driven  by  different  pump  configurations.  This  data  was 
collected  for  the  end  purpose  of  developing  effective  model  updating  strategies  for  the  detection  of  pipe 
rupture  using  networks  of  w  ireless  sensing  prototypes. 

5.2. 1.  Phase  II  Objectives: 

The  objectives  of  the  Phase  11  field  work  can  be  defined  as  follows: 

1.  Interface  Narada  wireless  sensors  with  the  flowmeters  already  installed  on  the  ONR  tabletop 
chilled  water  demonstration  system  using  suggested  improvements  from  Phase  1.  Also,  interface 
Narada  wireless  sensors  with  newr  flowmeters  installed  at  each  rupture  location. 

2.  Utilizing  9  different  pump  configurations,  collect  large  amounts  of  full-system  flow'  data  from  the 
tabletop  demonstrator  in  varying  states  of  damage  (undamaged,  rupture  1,  rupture  2,  rupture  1 
and  rupture  2). 

5.2.2.  Large  Scale  Data  Collection  for  Model  Updating  and  Plant  Damage  Detection: 

When  developing  the  model  updating  methodology  presented  in  Section  4,  it  was  found  that  in  order  to 
accurately  match  experimental  results  with  analytical  model  projections,  we  need  to  collect  data  over 
several  different  input-output  combinations.  Because  we  have  control  over  the  pumps  in  the  chilled  w^atcr 
environment,  it  is  decided  to  use  different  combinations  of  pumps  to  create  multiple  input-output 
scenarios.  Each  configuration  is  created  by  powering  dow  n  a  certain  subset  of  pumps.  The  configurations 
used  in  this  study  are  as  follows  (pump  locations  can  be  seen  in  Figure  5.1b): 

Configuration  #  Pumps 

1  1,2, 3,4 

2  1,2,3 

3  2,3,4 
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Ruptures  1  &  2  Rupture  2  Rupture  1  No  Rupture 


4  3,4 

5  1,3 

6  2,3 

7  1 

8  2 

9  3 

After  implementing  the  ehanges  suggested  in  Phase  1,  a  large  amount  of  clean  data  was  collected  from  the 
demonstrator.  A  typical  example  of  flowmeter  voltage  data  from  a  run  in  Pump  Configuration  #1  is 
shown  in  Figure  5.10.  When  this  data  is  converted  to  flowrate  using  the  technique  developed  in  Phase  1,  a 
rolling  average  of  flow  data  can  be  created.  Over  a  sufficiently  large  sample  of  data,  this  average  can  be 
used  as  a  representation  of  flow  at  that  point  in  time.  A  visual  example  of  the  Narada’s  final  flowrate 
output  (sampled  at  200  Hz)  can  be  seen  in  Figure  5.11.  The  blue  dots  in  this  figure  represent  the 
converted  and  sampled  flowrates  and  the  red  line  represents  the  rolling  average. 

Table  5.1  contains  the  numerical  representation  of  a  subset  the  averaged  flow  data  collected  in  Phase  11. 
For  eaeh  pump  configuration,  the  flows  recorded  at  each  of  the  ten  flowmeters  (FI  through  F10)  can  be 
seen.  Flow  meter  locations  can  be  found  in  Figure  5.1  (flowmeters  are  represented  by  pink  circles). 


Table  5.1.  Example  Phase  II  averaged  flowrate  data 


Pump  Confia 

FI 

F2 

F3 

£4 

£5 

F6 

F7 

F8 

F9 

F10 

Configuration  1 

13.37 

14.14 

21.2 

20.35 

1373 

1388 

18.67 

18.89 

12,41 

11.41 

Configuration  2 

10.53 

11.28 

24.23 

23.13 

10.5 

11.09 

21.68 

0 

9.54 

8.89 

Configuration  3 

9.8 

10.88 

0.58 

23.37 

10.5 

10.47 

21.66 

22.37 

9.63 

8.73 

Configuration  5 

6.94 

7.51 

27.59 

0 

7.09 

7.53 

24.43 

0 

6.58 

5.62 

Configuration  7 

4.06 

4.27 

30.34 

0 

3.75 

4.22 

0 

0 

3.46 

2.4 

Configuration  8 

3.41 

4.67 

0 

26.96 

3.67 

4.09 

0 

0 

3.36 

2.34 

Configuration  1 

13.48 

14.34 

21.52 

20.23 

13.85 

14.06 

18.93 

19.36 

12.45 

11.48 

Configuration  2 

10.61 

11.42 

24.52 

23.21 

10.59 

11.27 

22.05 

0.17 

9.51 

8.9 

Configuration  3 

9.94 

11.1 

0.53 

23.33 

10.75 

1081 

21.88 

22  75 

9.74 

8.6 

Configuration  5 

7.09 

7.55 

27.89 

0 

7.13 

7.66 

24.92 

0 

6.59 

5.74 

Configuration  7 

4.03 

4.18 

30.41 

0 

3.67 

4.11 

0 

0 

3.37 

2.38 

Configuration  8 

3.5 

4.78 

0 

26.83 

3.76 

4.16 

0 

0 

3.4 

2.45 

Configuration  1 

13.43 

14.39 

21.18 

19.89 

13.73 

13.92 

18.93 

19.57 

12.35 

11.33 

Configuration  2 

10.49 

11.34 

24.42 

22  93 

10.56 

11.06 

21.99 

0.2 

9.44 

8.74 

Configuration  3 

9.88 

11.03 

0.51 

23.07 

1063 

10.6 

21.8 

22.75 

9.63 

8.66 

Configuration  5 

7  01 

7.58 

27.66 

0 

7.31 

7.76 

24.82 

0.04 

6.73 

5.26 

Configuration  7 

4.12 

4.33 

30.24 

0 

3.84 

4.33 

0 

0 

3.55 

2.33 

Configuration  8 

3.33 

4.61 

0 

27.24 

3.65 

4.09 

0 

0 

3.36 

2.16 

Configuration  1 

13,48 

14.38 

21.55 

20.23 

13.78 

14.03 

19.04 

19.6 

12.63 

11.52 

Configuration  2 

10,48 

11.24 

24.79 

23.39 

10.51 

11.13 

22.09 

0.18 

9.62 

8.77 

Configuration  3 

9.97 

11.01 

0.48 

23.31 

10.69 

10.73 

21.92 

22.89 

9.88 

8.81 

Configuration  5 

7.11 

7.53 

27.85 

0 

7.22 

7.68 

25.02 

0 

6.83 

5.73 

Configuration  7 

4.06 

4.23 

30.33 

0 

3.66 

4.12 

0.85 

0 

3.46 

2.33 

Configuration  8 

3.33 

4.65 

0 

27.18 

3.59 

3.97 

0 

0 

3.39 

2.25 
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Experimental  Data  Run  6 
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Figure  5.10.  Typical  Phase  11  flowmeter  voltage  data,  sampled  at  10,000  Hz. 
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Figure  5.11.  Typical  Phase  II  converted  flowrate  data  (blue)  and  rolling  average  (red),  sampled  at 

200  Hz. 
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6.  Model  Updating  of  a  Chilled  Water  Plant  with  Rupture  Conditions: 

In  a  naval  environment,  the  ability  to  detect  distinct  changes  in  the  operational  condition  of  a  engineering 
plant  is  extremely  desirable.  In  this  study,  a  network  of  wireless  sensors  is  used  in  conjunction  with  an 
analytical  model  of  a  pipe  network  in  order  to  detect  damage  and/or  system  bifurcations  in  a  shipboard 
chilled  water  system.  To  this  end,  model  updating  methods  have  been  developed  and  embedded  within 
the  computational  core  of  the  Narada  wireless  sensing  unit.  These  methods  function  by  “updating" 
unknown  system  parameters  within  an  analytical  system  model  such  that  they  match  experimental  data 
collected  by  the  wireless  sensors.  Over  time,  if  and  when  these  updated  analytical  parameters  change, 
they  can  serve  as  indicators  of  the  severity  and  location  of  any  acquired  system  damage. 

For  the  purpose  of  this  study,  an  analytical  model  of  a  pipe  network  is  created  that  can  be  used  to 
determine  flowrate  information  in  each  pipe,  given  a  certain  pipe  and  pump  configuration.  This  analytical 
model  is  designed  using  the  linear  theory  method  as  proposed  by  Jeppson  (1976).  Essentially,  this  method 
is  used  to  satisfy  a  system  of  equations  where  an  independent  continuity  equation  (XQ,  =  0)  is  written  at 
each  pipe  junction  and  a  head  loss  equation  (XK;Q,  =  AH)  is  written  around  each  loop  in  the  network. 
Additionally,  equations  relating  pump  flow'  rate  and  pump  head  loss  arc  written  for  each  pump.  Using  an 
iterative  convergence  process,  a  set  of  system  flowrates  can  be  settled  upon  that  satisfy  head  loss,  pump, 
and  flowrate  continuity  requirements. 

To  validate  the  ability  of  this  model  updating  system  to  detect  and  locate  damage  in  a  chilled  water 
system,  the  ONR  tabletop  chilled  water  system  demonstrator  is  chosen  as  a  viable  testbed.  This  system 
contains  42  pipes,  4  pumps,  6  loads,  26  valves,  and  2  rupture  valves.  Each  of  the  10  pipes  containing 
either  a  load  or  a  pump  is  monitored  w'ith  a  flowmeter.  Because  of  the  ability  of  the  demonstrator  to 
simulate  pipe  rupture  through  the  use  of  two  rupture  valves,  pipe  rupture  is  chosen  as  the  primary  damage 
mechanism  for  investigation.  As  such,  a  properly  functioning  model  updating  system  will  be  capable  of 
detecting  when  one  (or  both)  of  the  rupture  valves  is  open,  and  howf  much  water  is  flowing  through  each 
of  them  given  a  set  of  flowmeter  data.  A  schematic  drawing  of  the  ONR  tabletop  demonstrator  can  be 
found  in  Figure  6.1 . 

6.1.  Examination  of  the  Proposed  Analytical  Model: 

In  the  initial  analytical  model  used  in  this  study  to  represent  the  flow  distribution  throughout  an 
interconnected  network  of  pipes,  the  Darcy-Wcisbach  equation  is  utilized  to  determine  head  loss  in  a  pipe 
with  a  given  length,  diameter,  flowrate,  and  friction  factor.  The  Darcy-Weisbach  equation  is  the  most 
accurate  known  method  for  modeling  flow  through  a  pipe,  but  it  relies  on  the  calculation  of  friction  factor 
from  a  set  of  equations  that  arc  dependent  on  the  type  of  flow  being  seen  in  the  pipe.  As  such,  these 
equations  are  highly  accurate  when  dealing  with  flow's  that  are  either  fully  laminar  or  fully  turbulent. 
Unfortunately,  because  of  the  scale  of  the  ONR  tabletop  chilled  water  system,  most  of  the  pipes  in  this 
network  fall  in  the  transition  range  between  laminar  and  turbulent  flow.  Therefore,  it  is  found  that  the 
Darcy-Wcisbach  equation  does  not  produce  consistent  results  across  the  board,  and  in  many  cases  it  may 
predict  flowrates  in  some  pipes  that  are  orders  of  magnitude  larger  than  what  is  seen  in  reality. 

As  a  result,  it  is  decided  to  utilize  the  Hazcn-Williams  equation  in  lieu  of  the  Darcy-Wcisbach  equation 
for  the  model  updating  work  performed  in  this  study.  The  Hazcn-Williams  equation  presents  a  much  more 
simplistic  (and  therefore  less  accurate)  representation  of  the  relationship  between  flowrate  and  head  loss 
in  a  given  pipe,  but  it  is  linear  throughout  the  transition  region  between  laminar  and  turbulent  flow,  and 
thus  allows  us  to  look  more  confidently  at  systems  where  some  pipes  are  laminar,  some  arc  turbulent,  and 
some  are  in-between.  In  the  model  updating  case,  where  we  are  much  more  interested  in  system  changes 
over  time  than  in  the  correlation  between  a  given  model  and  reality,  the  Hazen-Williams  equation  can  be 
an  excellent  choice. 
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6.2.  Model  Updating  using  Pipe  Diameters: 

In  the  initial  model  updating  method  employed  in  this  study,  pipe  damage  is  designed  to  be  associated 
direetly  with  pipe  diameter.  In  this  methodology,  the  diameter  of  eaeh  pipe  is  treated  as  an  unknown, 
constrained  by  a  range  defined  between  0%  and  200%  of  the  pipe’s  actual  diameter.  Then,  these 
unknowns  are  utilized  as  the  updating  parameters  in  the  model  updating  problem,  and  a  seareh  is 
launched  to  find  the  set  of  42  pipe  diameters  whieh  produec  a  set  of  analytical  system  flowrates  that 
match  most  elosely  the  10  experimental  flowrates  measured  in  the  tabletop  demonstrator.  The  thought 
behind  this  methodology  is  that,  in  the  undamaged  system,  the  set  of  diameters  that  the  model  updating 
algorithm  converges  upon  should  be  more  or  less  identical  to  the  aetua  pipe  diameters  in  the  physical 
system.  Then,  as  rupture  (or  bloekage)  is  introduced  in  a  given  pipe,  a  new  model  update  should  show  an 
increase  (or  decrease)  in  updated  pipe  diameter  for  that  given  pipe. 

It  is  quickly  found  that  by  using  this  methodology,  the  model  updating  algorithm  is  capable  of  finding  a 
set  of  pipe  diameters  that  produced  analytical  flowrates  exaetly  equivalent  to  those  measured 
experimentally.  Unfortunately,  the  updated  set  of  pipe  diameters  does  not  reflect  the  physical  system  in 
two  very  distinet  ways: 

1 .  There  is  no  correlation  found  between  updated  pipe  diameter  and  physical  pipe  diameter. 

2.  Sueeessive  model  updates,  while  always  effectively  finding  a  system  that  prediets  analytical 
flowrates  equivalent  to  experimental  flowrate  measurements,  returns  vastly  different  sets  of 
updated  pipe  diameters  eaeh  time  the  algorithm  is  run. 

As  such,  this  methodology  is  found  to  be  worthless  from  the  damage  detection  perspective,  as  there  (I)  is 
no  way  to  direetly  correlate  the  analytical  model  (pipe  diameters)  to  the  physical  system,  and  (2)  there  is  a 
seemingly  infinite  number  of  solutions  to  the  pipe  diameter  model  updating  problem. 

6.3.  Model  Updating  using  Hazen-Williams  Coefficients: 

In  an  attempt  to  improve  upon  the  initial  updating  methodology,  it  is  deeided  to  try  and  correlate  pipe 
damage  not  to  a  pipe’s  diameter,  but  to  a  pipe’s  Hazen-Williams  coefficient  (Chw)-  The  thought  behind 
this  methodology  is  that  beeause  the  Hazen-Williams  coefficient  is  tied  direetly  to  a  pipe’s  frietional 
properties  (and  thus  its  influence  on  head  loss),  it  may  better  refleet  system  ehanges  in  the  ease  of  a  pipe 
rupture  or  bloekage.  Using  this  methodology,  the  Hazen-Williams  eoeffleient  of  eaeh  pipe  is  treated  as  an 
unknown  model  updating  parameter,  constrained  by  a  range  defined  between  I  (very  rough  pipe)  and 


Figure  6.1.  Pipe  numbering  and  flow  direction  schematic  for  chilled  water  system  with  rupture 

valves  closed. 
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1000  (very  smooth  pipe).  Then,  a  seareh  is  launched  to  find  the  set  of  42  coefficients  which  produce  a  set 
of  analytical  system  flowrates  that  match  most  closely  the  10  experimental  flowrates  measured  in  the 
tabletop  demonstrator.  The  thought  behind  this  methodology  is  that,  in  the  undamaged  system,  the  set  of 
coefficients  that  the  model  updating  algorithm  converges  upon  should  correlate  to  the  actual  frictional 
capacity  of  each  pipe  in  the  physical  system.  Then,  as  a  rupture  (or  blockage)  is  introduced  in  a  given  pipe, 
a  newr  model  update  should  show  a  respective  increase  (or  decrease)  in  the  frictional  coefficient  for  that 
pipe. 

At  first  glance,  it  seems  as  if  this  methodology  is  preferable  to  the  updating  of  pipe  diameters,  as  Hazcn- 
Williams  coefficients  are  more  intuitively  tied  to  the  physical  actuality  of  friction  change  associated  w'ith 
pipe  rupture  or  pipe  blockage.  Unfortunately,  this  methodology  is  hampered  by  the  same  problems  found 
when  using  pipe  diameters  as  an  updating  parameter.  Specifically,  there  is  a  seemingly  infinite  number  of 
solutions  to  the  model  updating  problem,  making  damage  detection  very  difficult,  if  not  impossible. 

6.4.  Dimensionality  Reduction  of  the  Model  Updating  Problem: 

One  of  the  ways  in  which  the  problem  of  infinite  solutions  makes  itself  most  obvious  is  that  after  a 
completed  model  update,  two  pipes  which  have  the  same  length,  diameter,  and  internal  frictional  losses 
(bcnds/tees/valves)  may  be  found  to  have  vastly  different  Chw  values.  For  example,  pipes  9,  10,  41,  and 
42  (found  in  Figure  4)  should  have  very  similar  frictional  properties,  and  thus  very  similar  Chw  values. 
However,  it  is  found  solutions  obtained  using  the  model  updating  method  are  not  constrained  in  this  way. 
As  such,  it  is  decided  to  pursue  a  model  update  containing  only  12  different  Hazen-Williams  coefficients. 
One  coefficient  is  assigned  to  each  of  the  ten  pipes  containing  either  a  pump  or  a  load,  one  is  assigned  to 
each  of  the  four  pipes  with  tw  o  bends  (9,  10,  41,  and  42),  and  one  is  assigned  to  the  remainder  of  the 
pipes,  each  of  which  has  a  three-way  tee  at  each  end  and  no  bends. 

Immediately  after  applying  this  dimensionality  reduction  principle  to  the  model  updating  problem,  a  vast 
improvement  can  be  seen  in  the  quality  of  the  final  results.  Not  only  does  an  exact  match  between 
analytical  and  experimental  flowrates  continue  to  be  reached,  but  the  resulting  values  begin  to  show  a 
strong  correlation  with  what  would  be  expected  from  the  actual  physical  system.  For  example,  the  Chw 
value  in  the  smaller  pipes  (1,  3,  5,  7,  19,  28,  31,  34,  37,  40),  eaeh  of  w  hich  contain  several  obstructions,  is 
significantly  lowrer  (designating  more  friction)  than  the  Cfiw  values  found  for  the  larger,  unobstructed 
pipes.  Additionally,  a  marked  difference  is  seen  between  the  Chw  values  in  the  four  pipes  containing 
pumps  (1,  7,  31,  34)  than  in  those  with  loads.  Overall,  pipes  that  should  have  similar  properties  do  have 
similar  properties 

Unfortunately,  while  this  dimensionality  reduction  seems  to  improve  the  tie  between  the  physical  system 
and  the  model-updated  results,  variability  in  the  results  between  model  updating  runs  leads  to  the  belief 
that  the  problem  of  infinite  solutions  is  still  present  in  the  system. 

6.5.  Simultaneous  Updating  of  Multiple  Input-Output  Relations: 

Another  of  the  ways  in  which  the  problem  of  infinite  solutions  manifests  itself  in  the  damage  detection 
architecture  is  by  introducing  solutions  (i.c.  sets  of  Hazen-Williams  coefficients)  that  may  produce 
accurate  model  predictions  for  one  particular  input  (i.e.  all  four  pumps  engaged),  but  vastly  incorrect 
model  predictions  for  another  input  (i.e.  only  three  pumps  engaged).  As  such,  it  is  useful  to  use  multiple 
sets  of  flowmeter  data,  collected  using  different  input  parameters  (in  this  ease,  number  of  pumps  engaged) 
to  update  several  input-output  relations  simultaneously.  In  other  wwds,  by  utilizing  two  or  more  different 
sets  of  input-output  data  from  the  same  undamaged  system,  we  can  update  an  analytical  model  for  each 
type  of  input  using  one  unique  set  of  the  twelve  updating  parameters  discussed  in  the  previous  section. 
Hopefully,  this  method  will  decrease  significantly  the  number  of  possible  CHw  combinations  that  produce 
accurate  matches  between  analytical  and  experimental  results. 
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By  simultaneously  updating  multiple  analytical  models  of  different  input  scenarios,  it  is  found  that  while 
each  additional  model  slightly  degrades  the  accuracy  of  the  match  between  experimental  measurements 
and  analytical  predictions  (likely  due  to  experimental  error  or  model  imperfection),  the  updated  set  of 
Hazen-Williams  coefficients  becomes  largely  unique.  For  example,  when  the  sets  of  data  generated  with 
three  pumps  and  four  pumps  are  updated  simultaneously,  a  unique  set  of  10  of  the  12  updated  Hazen- 
Williams  coefficients  can  be  found  (those  corresponding  to  the  monitored  pipes)  that  predicts  flowrates  to 
w  ithin  0.01  mL/see  of  all  recorded  measurements.  However,  the  two  Hazen-Williams  coefficients  that  are 
applied  to  the  remainder  of  the  pipes  in  the  system  continue  to  fluctuate,  leaving  a  still  infinite  number  of 
possible  system  configurations. 

6.6.  Addition  of  One  Flowmeter  to  Ensure  Asymmetric  Monitoring: 

Upon  examining  why  only  the  Hazen-Williams  coefficients  associated  with  monitored  pipes  could  be 
determined  using  the  previous  approach  with  multiple  input-output  relations,  it  became  apparent  that 
knowing  full  information  about  these  ten  pipes  did  not  necessarily  imply  anything  about  the  remainder  of 
the  pipe  network.  In  other  words,  even  if  flowr  rates  and  head  losses  in  these  monitored  pipes  are  fixed,  the 
flow  rates  and  head  losses  in  the  surrounding  pipes  can  fluctuate  at  will,  due  to  the  symmetrical  nature  of 
the  pipe  network  and  the  centralized  positioning  of  the  existing  flow  meters.  As  such,  it  is  proposed  that  an 
additional  flowmeter  be  installed  in  one  of  main  pipes  used  to  connect  parts  of  the  system  (i.e.  pipe  20,  21, 
22,  or  23).  Potentially,  this  w  ill  force  the  remainder  of  the  system  into  one  unique  flow  configuration,  and 
as  such,  we  should  be  able  to  detect  one  unique  set  of  Hazen-Williams  coefficients. 

After  updating  the  same  set  of  models  used  in  the  previous  section,  it  is  found  that  the  addition  of  one 
more  flowmeter,  installed  in  pipe  21,  does  in  fact  allow  the  model  updating  process  to  detect  one  unique 
set  of  Hazen-Williams  coefficients  for  a  given  system  state.  The  unique  nature  of  this  solution  is 
conducive  to  damage  detection  applications,  where  a  solid  baseline  value  is  necessary  in  order  to 
effectively  look  for  changes  in  system  properties. 

6.7.  Model  Updating  for  Pipe  Rupture  Detection  and  Quantification: 

Having  discovered  a  method  of  effectively  determining  system-wide  properties  through  model  updating 
(by  adding  an  additional  flowmeter,  using  dimensionality  reduction  and  by  simultaneous  updating 
multiple  input-output  scenarios),  it  is  now  possible  to  look  into  the  detection  and  location  of  damage  (pipe 
rupture)  within  the  tabletop  demonstrator.  In  order  to  most  accurately  model  the  experimental  situation 
present  in  the  physical  tabletop  demonstrator,  it  is  decided  to  model  the  ruptures  exactly  as  they  exist  in 
reality.  In  other  words,  a  pipe  is  added  at  the  location  of  each  rupture  valve  which  joins  the  main 
demonstrator  pipe  network  to  the  reservoir  which  provides  water  to  the  system.  Thus,  in  order  to 
accommodate  the  rupture  valves  and  the  reservoir,  the  number  of  pipes  is  increased  from  42  to  47.  This 
new  analytical  model  can  be  seen  in  Figure  6.2.  Ideally,  an  additional  flowmeter  would  be  placed  in  pipe 
21 ,  but  due  to  demonstrator  restrictions,  this  addition  was  physically  impossible. 

Given  the  analytical  model  proposed  above,  it  is  postulated  that  rupture  can  be  detected  and  localized 
(given  a  previous  know  ledge  of  the  two  rupture  locations)  by  using  the  follow  ing  procedure: 

1.  Collect  flowmeter  data  in  the  undamaged  (i.e.  rupture  valves  closed)  case  for  two  or  more 
different  input  scenarios  (created  by  varying  the  pump  configuration). 

2.  Simultaneously  update  all  12  unknown  Chw  values  using  the  analytical  models  corresponding  to 
the  varying  pump  configurations,  while  forcing  Chw  in  pipes  44  and  46  to  be  equal  to  0.01 .  (This 
simulates  near-infinite  friction,  and  thus  a  flow  situation  where  both  rupture  valves  are  closed.  A 
Chw  of  0  cannot  be  used  because  of  divide-by-zero  issues.)  The  resulting  updated  Chw  values 
represent  the  baseline,  undamaged  state  of  the  demonstrator. 

3.  Collect  flowmeter  data  in  all  three  undamaged  states  (i.e.  one  or  both  rupture  valves  open)  for 
each  of  the  pump  configurations. 
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4.  For  each  damage  state,  simultaneously  update  eaeh  of  the  appropriate  analytical  models 
(corresponding  to  pump  configuration)  using  the  CH w  value  in  pipes  44  and  46  as  the  updating 
parameters.  The  baseline,  undamaged  Chw  values  found  in  step  2  are  utilized  for  all  other  pipes. 
The  resulting  updated  Chw  values  for  pipes  44  and  46  represent  the  ruptured  state  of  the 
demonstrator.  The  larger  these  values  are  relative  to  the  baseline  ease  (Step  2),  the  greater  the 
degree  of  rupture  in  the  associated  location. 

6.8  Experimental  Model  Updating  Results  for  Pipe  Rupture  Deteetion  and  Quantification: 

As  described  in  Seetion  5,  during  Phase  II  experimental  testing  a  large  set  of  data  is  collected  using  a 
wide  range  of  pump  configurations  to  drive  the  system.  While  all  nine  pump  configurations  were  tested,  it 
was  found  that  four  of  the  configurations  were  better  suited  to  driving  water  into  the  open  rupture  valves. 
As  sueh,  it  was  deeided  to  only  use  these  four  pump  configurations  (1,  2,  3,  and  5)  for  the  purpose  of 
model  updating-based  damage  deteetion.  In  eaeh  of  these  pump  configurations,  data  was  eolleeted  for  the 
undamaged  system,  the  system  with  rupture  valve  1  open,  the  system  with  rupture  valve  2  open,  and  the 
system  with  both  rupture  valves  open. 

Before  any  experimental  data  is  used  to  test  the  damage  deteetion  strategy  proposed  herein,  it  is  deeided 
to  utilize  a  set  of  simulated  data  derived  from  the  analytical  model  seen  in  Figure  6.2  and  explained  above. 
As  sueh,  the  model  is  driven  with  the  same  set  of  four  pump  configurations,  and  simulated  data  is 
eolleeted  in  the  undamaged  system,  the  system  with  rupture  valve  1  open,  the  system  with  rupture  valve  2 
open,  and  the  system  with  both  rupture  valves  open.  While  several  different  objective  functions  were 
effectively  utilized,  it  was  found  that  the  most  effective  objective  function  consisted  of  the  average 
analytical  difference  between  analytical  model  flowrates  and  experimental  flowmeter  data  (essentially 
Equation  4.3). 

Using  the  four  different  input  eases  generated  by  varying  the  pump  configuration,  it  is  found  that  the 
proposed  model  updating  procedure  is  capable  of  not  only  detecting,  but  also  quantifying  (to  an  extent) 
rupture  in  the  two  locations  targeted  in  this  study.  Cases  are  tested  where  rupture  valve  1  is  open,  where 
rupture  valve  2  is  open,  and  where  both  rupture  valves  are  open.  In  all  eases,  the  model  updating 
algorithm  is  capable  of  detecting  damage.  However,  the  degree  to  which  that  damage  can  be  loeated  and 
quantified  varied  from  ease  to  ease. 


Figure  6.2.  Pipe  numbering  and  flow  direction  schematic  for  chilled  water  system  with  rupture 

valves  open. 
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Beginning  with  the  simulated  set  of  data,  the  12  unknown  CHw  values  are  updated  in  an  effort  to  find  an 
analytical  system  that  produces  flowrate  data  that  is,  on  average,  as  elose  as  possible  to  the  simulated  data 
from  the  undamaged  chilled  water  system.  The  CHw  values  for  pipes  44  and  46  are  set  to  0.01  as 
explained  above.  As  seen  in  Figure  6.3,  the  model  updating  procedure  is  able  to  find  a  set  of  CHw  values 
that  ereate  analytical  flowrates  very  elose  to  those  produced  in  the  (undamaged)  simulated  system  (within 
1%  error).  These  baseline  values  are  then  used  to  update  the  possible  damage  locations  (rupture  1  and 
rupture  2)  in  each  of  the  four  possible  damage  scenarios.  It  ean  be  seen  in  Figure  6.4,  Figure  6.5,  Figure 
6.6,  and  Figure  6.7  that  CHw  values  for  pipes  44  and  46  ean  be  found  that  produce  analytical  flowrates 
very  elose  to  those  found  in  the  corresponding  simulated  systems  (within  2%  error  in  all  eases).  When 
looking  at  how  the  updated  Chw  values  for  pipes  44  and  46  relate  to  their  respective  damage  eases  (Figure 
6.8),  it  ean  be  elearly  seen  that  the  damaged  scenarios  (where  one  or  both  of  the  rupture  valves  are  open) 
show  very  distinctive  signs  of  damage  when  compared  with  the  undamaged  scenario.  However,  while 
there  does  seem  to  be  a  correlation  between  whieh  the  physical  and  the  updated  rupture  locations  (Chw 
values  arc  higher  for  the  physical  rupture  locations),  the  faet  that  the  model  updating  procedure  diagnoses 
some  degree  of  rupture  in  both  pipes  for  all  damage  scenarios  (even  when  only  one  rupture  location 
exists),  this  method  appears  questionable  for  locating  and  quantifying  damage  within  a  pipe  network. 

Moving  to  the  set  of  experimental  data,  the  12  unknown  Chw  values  are  again  updated  in  an  effort  to  find 
an  analytical  system  that  produces  flowrate  data  that  is,  on  average,  as  elose  as  possible  to  the 
experimental  data  from  the  undamaged  chilled  water  system.  The  Chw  values  for  pipes  44  and  46  arc  set 
to  0,01  as  before.  As  seen  in  Figure  6.9,  the  model  updating  procedure  is  able  to  find  a  set  of  Chw  values 
that  ereate  analytical  flowrates  reasonably  close  to  those  produced  in  the  (undamaged)  experimental 
system  (within  20%  error  in  all  eases  and  mueh  eloser  in  most  eases).  These  baseline  values  are  then  used 
to  update  the  possible  damage  locations  (rupture  1  and  rupture  2)  in  eaeh  of  the  four  possible  damage 
scenarios.  It  can  be  seen  in  Figure  6.10,  Figure  6.11,  Figure  6.12,  and  Figure  6.13  that  Chw  values  for 
pipes  44  and  46  can  be  found  that  produce  analytical  flowrates  very  elose  to  those  found  in  the 
corresponding  simulated  systems  (again,  within  20%  error  in  all  eases  and  much  closer  in  most  cases). 
When  looking  at  how  the  updated  Chw  values  for  pipes  44  and  46  relate  to  their  respective  damage  cases 
(Figure  6.14),  it  can  be  seen  (although  somewhat  less  elearly  than  in  the  case  of  simulated  data),that  the 
damaged  scenarios  (where  one  or  both  of  the  rupture  valves  are  open)  show  some  signs  of  damage  w  hen 
compared  with  the  undamaged  scenario.  There  even  seems  to  be  some  degree  of  correlation  between 
which  the  physical  and  the  updated  rupture  locations  (Chw  values  are  higher  for  the  physical  rupture 
locations),  especially  in  the  ease  of  rupture  1  and  ruptures  1  &  2.  However,  the  faet  that  the  model 
updating  procedure  does  not  diagnose  a  greater  degree  of  rupture  at  rupture  location  2  for  the  case  of 
rupture  2,  this  method  (as  with  the  simulated  data)  appears  questionable  for  locating  and  quantifying 
damage  within  a  pipe  network. 
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Figure  6.3.  Model  updating  accuracy  for  a  baseline  update  using  simulated  data. 
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Figure  6.4.  Model  updating  accuracy  for  an  unruptured  update  using  simulated  data 
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Figure  6.5.  Model  updating  accuracy  for  a  rupture  case  1  update  using  simulated  data. 
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Figure  6.6.  Model  updating  accuracy  for  a  rupture  case  2  update  using  simulated  data. 
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Figure  6.7.  Model  updating  accuracy  for  a  rupture  ease  1&2  update  using  simulated  data. 
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Figure  6.8.  Model  updating  damage  detection  results  for  all  damage  cases  using  simulated  data. 
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Updated  (Baseline)  vs.  Experimental  Flowrates 
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Figure  6.9.  Model  updating  accuracy  for  a  baseline  update  using  experimental  data. 
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Figure  6.10.  Model  updating  accuracy  for  an  unruptured  update  using  experimental  data. 
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Figure  6.11.  Model  updating  accuracy  for  a  rupture  case  1  update  using  experimental  data. 
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Figure  6.12.  Model  updating  accuracy  for  a  rupture  case  2  update  using  experimental  data. 
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Damage  Case 


Updated  (Rupture  1&2)  vs.  Experiment®*  Flowrates 
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Figure  6.13.  Model  updating  accuracy  for  a  rupture  case  1&2  update  using  experimental  data. 


Experimental  Model  Updating  Results 


Figure  6.14.  Model  updating  damage  detection  results  for  all  damage  cases  using  experimental  data. 
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7.  Market-based  computational  task  assignment  within  a  WSN: 

Having  developed  a  framework  in  which  a  wireless  sensing  network  can  execute  complex  engineering 
analyses  (such  as  model  updating)  in  a  distributed  fashion,  the  next  logical  step  is  to  expand  this 
computing  paradigm  to  allow  a  sensor  network  to  perform  multiple  computational  tasks  simultaneously. 
However,  one  of  the  key  challenges  inherent  in  the  simultaneous  deployment  of  multiple  distributed 
algorithms  on  a  wireless  sensing  network  is  that  in  the  wireless  environment,  many  of  the  system 
resources  (such  as  battery  power,  data  storage  capacity,  MPU  time,  wireless  bandwidth,  etc.)  required  to 
perform  complex  computational  tasks  are  available  only  in  a  limited  manner.  As  such,  it  is  important  to 
devise  a  method  of  optimally  distributing  and  consuming  these  scarce  system  resources  throughout  the 
network.  One  approach  to  this  problem  is  to  apply  free-market  economics  to  help  allocate  these  resources. 
Free-market  economies  can  be  thought  of  as  large  collections  of  autonomous  market  agents  (participants) 
such  as  producers  (sellers)  and  consumers  (buyers),  among  others.  In  this  free-market  framework,  each 
agent  is  forced  to  compete  against  other  agents  in  a  competitive  marketplace  with  scarce  resources. 

In  the  case  of  distributing  computational  tasks  throughout  wireless  sensing  networks,  a  market  can  be 
formed  where  the  only  service  for  sale  is  a  block  of  CPU  time  w  hich  can  be  devoted  to  any  computational 
task.  For  our  purposes,  we  will  deal  with  a  network  in  which  each  computational  task  is  a  simulated 
annealing  search  problem  of  complexity  Csa  that  can  be  broken  into  separate  blocks  of  N$a  SA  search 
iterations.  At  any  given  instant  in  time,  a  certain  number  of  processors,  PSa ,  arc  already  associated  with 
each  task,  and  each  task  has  already  seen  search  progress  made  at  SA  temperature  steps  up  to  and 
including  TSSa-  In  order  to  provide  a  viable,  robust  testbed  in  which  to  work  out  and  validate  the  market- 
based  concepts  proposed  herein,  the  w-Quccns  problem  is  chosen  as  a  simple,  easy  to  implement 
combinatorial  optimization  problem  that  will  fit  nicely  within  the  distributed  model  updating  framework 
that  was  developed  in  Section  4  and  applied  to  a  chilled  water  system  in  Section  6. 

7.1  The  n-Queens  problem: 

The  tt-Quecns  problem  is  a  well-known  benchmark  problem  for  evaluating  the  performance  (i.e.  speed 
and  efficiency)  of  optimization  or  search  algorithms  developed  in  the  computer  science  community.  The 
objective  of  the  n-Qucens  problem  is  to  place  n  chess  queens  on  an  n  x  n  chess  board  (where  n  >  4)  such 
that  no  queen  can  attack  another  queen  following  basic  chess  rules.  In  other  words,  no  queen  can  be 
placed  on  the  same  row,  column,  or  diagonal  as  another  queen.  One  of  many  solutions  to  the  8-Quccns 
problem  can  be  seen  in  Figure  7.1(b). 

The  /7-Qucens  optimization  problem  proceeds  by  attempting  to  minimize  an  objective  function,  £,  which 
sums  the  number  of  conflicts  betw  een  queens  in  a  given  chess  board  configuration.  In  an  analytical  sense, 
if  a  queen  is  at  position  (/,  J),  it  is  in  direct  conflict  with  any  queen  at  position  (/,/),  where  i  =  I  (same 
column),  or  j  =  J  (same  row),  or  |  i  —  I\=\j  -J\  (same  diagonal).  So,  if  we  let  qtj  represent  each  square  on  a 
chess  board,  and  if  wc  set  qy  equal  to  1  if  there  is  a  queen  at  position  (/,y),  and  0  otherw  ise,  we  can  create 
an  appropriate  objective  function,  £,  as  follows: 
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Figure  7.1.  (a)  Initial  board  configuration  for  8-Queens  problem  and 
(b)  optimal  solution  for  the  8-Qucens  problem 
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with  the  first  term  summing  row  conflicts,  the  second  term  summing  column  conflicts,  the  third  term 
summing  upper  diagonal  conflicts,  and  the  fourth  term  summing  lower  diagonal  conflicts.  Each 
combination  of  squares  qy  and  qu  returns  1  if  there  is  a  queen  conflict  and  0  if  there  is  not,  leading  to  a 
sum  equal  to  the  total  number  of  conflicts.  To  minimize  duplicate  conflicts,  each  square  on  the  chess 
board  is  evaluated  only  once  against  all  other  squares. 

For  the  implementation  of  the  //-Queens  problem  in  this  study,  wc  choose  to  start  with  a  board 
configuration  such  that  a  queen  is  placed  on  each  diagonal  square  (/,/),  i.e.  where  /  =  j,  as  seen  in  Figure 
7.1(a).  Clearly,  in  this  initial  state,  each  queen  is  in  conflict  with  all  other  queens.  New  search  states  arc 
generated  by  swapping  the  queens  laying  on  two  randomly  selected  rows,  while  retaining  each  queen’s 
initial  column.  In  this  way,  there  is  always  one  queen  in  each  row  and  one  queen  in  each  column.  This 
state  generation  method  allows  for  significantly  faster  convergence  of  the  optimization  problem,  as  the 
first  two  terms  of  the  objective  function  (Equation  7.1)  can  be  ignored.  The  //-Queens  problem  is  an  ideal 
testbed  in  this  study  because  it  allows  us  to  easily  explore  computational  tasks  of  varying  complexity  by 
simply  increasing  the  /7-Quccns  problem  size. 

7.2  Market-based  task  assignment: 

One  incredibly  complex  system  that  is  optimally  controlled  in  a  decentralized  manner  is  that  of  the  frcc- 
market  economy,  where  scarce  societal  resources  are  distributed  based  on  the  local  interactions  of  buyers 
and  sellers  who  obey  the  laws  of  supply  and  demand.  Recently,  researchers  have  begun  to  utilize  market- 
based  concepts  for  the  control  or  optimization  of  complex  systems,  most  often  in  the  realm  of  computer 
architecture  where  a  market  analogy  is  useful  for  modeling  computer  systems  such  as  memory  usage  or 
network  traffic.  Perhaps  the  greatest  benefit  of  market-based  optimization  is  that  can  often  yield  a  Pareto 
optimal  solution.  Pareto  optimal  can  be  defined  by  a  market  in  a  competitive  equilibrium  where  no  market 
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Figure  7.2.  Buyer/sellcr  distinction  in  the  market-based  task  assignment  model 


participant  can  reap  the  benefits  of  higher  utility  or  profits  without  causing  harm  to  other  participants 
when  a  resource  allocation  change  is  made. 

In  contrast  to  a  simple  auction-based  system  which  could  be  used  to  crudely  optimize  the  market  based  on 
only  one  goal,  such  as  minimizing  time  to  completion,  a  more  robust  market-based  task  assignment 
scheme  can  be  developed  that  can  optimally  allocate  resources  in  the  midst  of  several  market  goals.  This 
is  accomplished  through  the  use  of  buyer  and  seller  “utilities’'.  By  embedding  within  each  market  agent 
the  desire  to  optimize  an  individual  utility  function,  competing  goals  can  be  settled  through  market  means 
(supply,  demand,  price,  etc).  The  result  is  a  Pareto-optimal  allocation  of  scare  system  resources.  In  this 
example,  we  are  interested  in  three  distinct  (but  possibly  competing)  performance  objectives:  (1) 
completing  all  required  computational  tasks  as  quickly  as  possible,  (2)  minimizing  power  consumed  by 
the  sensor  network,  and  (3)  functioning  as  robustly  and  as  reliably  as  possible.  In  order  to  measure  the 
ability  of  the  market-based  system  to  address  these  three  competing  performance  objectives,  several 
performance  metrics  can  be  created  and  utilized:  (1)  the  time  required  to  complete  eaeh  task,  (2)  the 
number  of  wireless  transmissions  required  to  complete  each  task,  and  (3)  the  number  of  sensor  and 
communication  failures  encountered  during  each  task. 

7.2.1  B  uyer/sel  l  er  framework: 

As  seen  in  Figure  7.2,  the  sellers  in  this  market  can  be  defined  as  the  set  of  sensors  in  the  wireless 
network  not  currently  working  on  any  computational  task.  These  sensing  units  will  be  “selling”  their 
computational  abilities  to  a  number  of  buyers,  represented  by  the  set  of  sensors  most  recently  added  to 
eaeh  existing  computational  task.  For  example,  in  Figure  7.2  there  are  four  tasks,  so  there  are  four  buyers. 
In  order  to  address  all  three  goals  of  this  market  in  a  streamlined  manner,  buyers  and  sellers  arc  assigned 
different  objectives.  In  this  market,  all  buyers  work  both  to  minimize  the  overall  time  spent  computing 
and  to  maximize  network  reliability  by  (1)  minimizing  CPU  time  required  to  complete  each  task,  (2) 
minimizing  CPU  time  lost  due  to  sensor  failure,  and  (3)  minimizing  CPU  time  lost  due  to  communication 
failure.  Sellers,  on  the  other  hand,  work  to  minimize  network  power  consumption.  Beeause  the  wireless 
radio  consumes  significantly  more  power  than  any  other  sensor  component,  sellers  aim  to  minimize 
power  consumption  by  minimizing  the  number  of  wireless  communications  required  to  complete  eaeh 
task. 

7.2.2  Utility  function  development : 

In  light  of  this  framework,  it  is  necessary  to  first  outline  the  utility  functions  associated  with  both  buyers 
and  sellers.  These  utility  functions  will  govern  whether  or  not  the  buyer  for  each  computational  task  will 
plaee  a  bid  on  the  services  of  a  given  seller  and  which  buyer,  if  any,  a  seller  will  sell  its  services  to.  On 
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Figure  7.3.  (a)  Experimental  time  data  for  each  temperature  step/processor  combination,  (b) 
graphical  representation  of  experimentally  determined  T speed,  and  (c)  the  accompanying  analytical 

model  Fit 


the  buyer  side,  a  utility  function,  Ubuyer  can  be  intuitively  thought  of  as  the  total  amount  of  computational 
savings  a  computational  task  gains  by  adding  an  additional  processing  node,  and  can  be  defined  as 
follows: 


U  BW'ER  ~  ? SPEED  ~  a  BUYER  *  ^ SENSOR  FAIL  ~  ft  BUYER  ’  ? COMM  FAIL  {12) 

where  Tspeed,  Tsensorfail,  TCommfail ,  o.  buyer,  and  ft  buyer  are  defined  in  detail  below. 

For  any  simulated  annealing  task,  Tspeed  represents  the  decrease  in  expected  time  required  to  complete  the 
task  brought  about  by  the  addition  of  one  processor.  While  there  is  no  way  of  directly  formulating  an 
analytical  expression  for  this  value,  a  trend  can  be  established  for  any  computational  task  by  looking  at 
the  average  amount  of  time  it  takes  a  task  of  complexity  Csa  to  converge  on  a  solution  from  a  given 
temperature  step,  TSsa  with  a  given  number  of  processors,  Psa-  Using  data  collected  over  a  large  number 
of  experimental  trials,  the  time  saved  by  adding  an  additional  processor  at  any  given  search  state  can  be 
expressed  as  the  difference  between  the  average  time  required  to  complete  a  search  where  PSA  nodes  are 
currently  searching  up  to  temperature  step  TSsa  and  the  average  time  required  to  complete  a  search  where 
Psa+\  nodes  are  currently  searching  up  to  temperature  step  75&*+l .  Figure  7.3(a)  shows  the  time  required 
to  complete  the  100-Queens  problem  from  a  given  state  {Psa,  TSsa )  in  the  search,  and  Figure  7.3(b)  show's 
the  amount  of  time  saved  by  adding  an  additional  processing  node  w  hen  the  100-Queens  problem  is  in  a 
given  state  {Psa,  TSsa)-  Note  that  if  all  impossible  states  (w  here  PSA  >  TSsa)  are  excluded,  this  relationship 
is  independent  of  Psa,  and  can  be  approximated  by  an  easily  computable  algebraic  function: 

TSPEEo(CSAJSSAha  +  - . P^“-{p-y).e-TS‘* 

l  +  c  ;  (73) 

where  the  values  for  a,  /?,  S ,  and  y  are  specific  to  each  task  complexity,  Csa,  and  are  tabulated  in  Table  7.1. 
The  quality  of  the  analytical  fit  provided  from  this  function  for  the  1 00-Queens  problem  can  be  seen  in 
Figure  7.3(c).  Fits  of  similar  quality  can  be  found  for  all  other  problem  complexities  considered  (25- 
Quccns,  50-Queens,  and  75-Queens). 

For  any  simulated  annealing  task,  Tsensorfail  represents  the  increase  in  expected  processing  time  lost  due 
to  sensor  failure  brought  about  by  the  addition  of  one  processor.  Unlike  Tspeed ,  this  quantity  can  be 
derived  analytically.  Intuitively,  if  any  sensor  succumbs  to  either  hardware  or  software  failure  while  it  is 
involved  in  a  parallel  SA  task,  all  work  done  by  the  failed  node,  as  well  as  all  nodes  below  it  would  be 
lost.  As  such,  Tsensorfail  can  be  expressed  as  the  amount  of  time  required  for  the  newly  added  processor 
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to  complete  NSa  search  iterations  multiplied  by  the  probability  that  either  it  or  any  one  of  the  PSa 
processors  above  it  in  the  search  tree  succumbs  to  sensor  failure.  Analytically,  this  value  can  be  expressed 
as: 


1  SENSOR  FAIL 


(Psa)=(i- 


P  SENSOR  SUCCESS 


1 


(7.4) 


where  T^Sa  is  the  average  time  required  for  one  sensor  to  complete  NSa  search  iterations  and  Psensor 
success  is  the  probability  that  a  given  sensor  completes  its  NSa  search  iterations  without  failing.  This  value 
is  dependent  on  the  wireless  sensor  platform  being  used  and  the  environment  in  which  it  is  deployed,  but 
is  typically  very  high  (>0.999). 


For  any  simulated  annealing  task,  Tcommfail  represents  the  increase  in  expected  processing  time  lost  due 
to  communications  failure  brought  about  by  the  addition  of  one  processor.  Like  Tsensor  fail,  this  quantity 
can  also  be  derived  analytically.  In  a  similar  fashion,  if  any  sensor  loses  communication  with  its  parent 
while  it  is  involved  in  a  parallel  SA  task,  any  work  done  by  the  failed  node  and  all  nodes  below  it  would 
be  lost.  As  such,  Tcomm  fail  can  be  expressed  as  the  amount  of  time  required  for  the  newly  added 
processor  to  complete  N$a  search  iterations  multiplied  by  the  probability  that  either  it  or  any  one  of  the 
Psa-  1  processors  immediately  above  it  in  the  search  tree  loses  parental  communication.  At  first,  it  appears 
that  this  value  can  be  calculated  in  a  manner  similar  to  that  of  Tsensor  fail*  but  the  probability  of  success  in 
each  parent-child  communication  link  is  not  only  dependent  on  Psa,  but  also  on  the  signal  strength  of  the 
respective  w  ireless  communication  link.  This  value  can  be  quantified  for  each  communication  link,  c,  and 
represented  by  a  radio  signal  strength  indicator  (RSSIc)-  As  such,  an  analytical  value  for  Tcommfail  can  be 
expressed  as: 

Tcomm  fail  (Psa  >  KSSIC )  =  Pcomm  fail  '  ^nsa  (75) 

w'here  T^sa  is  as  before  and 

p«4 

_  ,  pi  Pcomm  success 
Pcomm  fail  -  1  ~  ]  [  .  -oa{4o.o+rssic) 

e  (7.6) 

Having  examined  in  more  detail  the  derivation  of  Tspeed,  TSensor  fail,  and  Tcommfail ,  it  can  now'  be  seen 
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from  Equation  7.2  that  aBVyER ,  and  Pbuyer  are  weighting  parameters  that  allow  the  wireless  sensor  network 
to  exactly  prioritize  between  speedup,  communication  reliability,  and  sensor  reliability.  This  type  of 
weighting  creates  an  extremely  adaptable  network  that  can  change,  in  real-time,  to  shifting  computing 
needs  within  a  wireless  network. 

On  the  seller  side,  a  somewhat  simpler  utility  function,  USeller ,  can  be  developed  in  a  similar  fashion  to 
U buyer-  Intuitively,  seller  utility  can  be  thought  of  as  the  total  amount  of  additional  power  a  computational 
task  requires  as  a  result  of  adding  an  additional  processing  node.  Because  the  majority  of  power 
consumption  in  a  wireless  sensing  device  comes  from  the  wireless  radio  (which  consumes  five  times 
more  power  than  a  microcontroller),  the  seller  can  maximize  its  utility  by  minimizing  the  amount  of  time 
the  wireless  network  spends  communicating.  As  such,  Useller  can  be  defined  as  follows: 

U  seller  =  ~~TComm  (7.7) 

For  any  simulated  annealing  task,  Tcomi  represents  the  increase  in  expected  communication  time  required 
to  complete  the  task  brought  about  by  the  addition  of  one  processor.  Much  like  TSpeed,  there  is  no  way  of 
directly  formulating  an  analytical  expression  for  this  value.  As  such,  a  communications  trend  can  be 
established  for  any  computational  task  by  looking  at  the  average  amount  of  time  that  a  task  of  complexity 
CSa  spends  communicating  before  converging  on  a  solution  from  a  given  temperature  step,  TSsa,  with  a 
given  number  of  processors,  PSa •  Using  data  collected  over  a  large  number  of  experimental  trials,  this 
value  can  be  expressed  as  the  difference  between  the  average  number  of  communications  required  to 
complete  a  search  wrhere  Psa  nodes  are  currently  searching  up  to  temperature  step  TSsa  and  the  average 
number  of  communications  required  to  complete  a  search  where  Psa+  1  nodes  are  currently  searching  up 
to  temperature  step  75^+1.  Figure  7.4(a)  shows  the  average  number  of  communications  required  to 
complete  the  100-Queens  problem  from  a  given  state  (Psa,  TSsa)  in  the  search,  and  Figure  7.4(b)  shows 
the  amount  of  additional  communications  required  by  adding  one  more  processing  node  w'hen  the  100- 
Quecns  problem  is  in  a  given  state  (Psa,  TSsa)-  Note  that  all  surface  plots  in  Figure  7.4  are  presented  on  a 
logarithmic  scale.  This  is  because  the  number  of  communications  required  during  the  first  temperature 
step  are  magnitudes  larger  than  what  is  required  during  any  other  temperature  step  and  a  log  scale  is 
required  in  order  to  properly  view  the  remainder  of  the  graph.  This  phenomenon  occurs  because  of  the 
large  number  of  new  global  minimum  states  that  are  found  at  the  beginning  of  the  SA  search  process. 
Using  simple  linear  approximations,  the  surface  in  Figure  7.4(b)  can  be  approximated  by  the  following 
easily  computable  algebraic  function: 


Figure  7.4.  (a)  Experimental  time  data  for  each  temperature  stcp/processor  combination,  (b) 
graphical  representation  of  experimentally  determined  Tcomi ,  and  (c)  the  accompanying  analytical 

model  fit. 
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and  where  values  for  a,  /?,  and  r5  are  specific  to  each  task  complexity,  Csa ,  and  arc  tabulated  in  Table  7.1. 
The  quality  of  the  analytical  fit  provided  by  this  function  for  the  100-Queens  problem  can  be  seen  in 
Figure  7.4(c).  Fits  of  similar  quality  can  be  found  for  all  other  problem  complexities  (25-Queens,  50- 
Queens,  and  75-Queens). 


Interestingly,  the  experimentally  derived  coefficients  a,  p ,  and  S  are  linearly  correlated  to  the  problem 
complexity,  CSa,  allowing  us  to  easily  expand  the  algebraic  expression  for  Tcomm^o  /7-Queens  problems  of 
any  arbitrary  non-negative  complexity: 
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7.2.3  W ire  less  tas  k  Jis  trihut  io t ;  algorithm : 

Having  developed  utility  functions  associated  with  both  buyers  and  sellers,  it  is  now  possible  to  create  a 
methodology  with  which  wireless  sensing  units  can  buy  and  sell  processing  time.  By  expanding  on  the 
fundamental  principles  of  an  auction,  the  following  procedure  is  developed:  (Note  that  at  any  time,  a  new 
computational  task  can  be  assigned  to  any  wireless  sensing  unit  not  otherwise  oeeupied). 

Step  1)  All  sensing  units  not  currently  processing  will  broadcast  their  availability  to  the 
network. 

Step  2)  Buyers  at  the  bottom  of  each  existing  computational  task  chain  will  calculate  U buyer 
based  on  the  computational  task  they  are  working  on,  and  place  a  bid  of  Ubuyer  if 
U buyer  >  0. 

Step  3)  Sellers  will  calculate  U seller  based  on  each  proposed  computational  job. 

Step  4)  Once  all  bids  have  been  received,  sellers  will  calculate  their  expected  profit  from 

each  proposed  job  using  a  market  power  /  speed  exchange  rate  iy market)  that 
represents  the  minimum  number  of  seconds  of  computational  speedup  that  must  be 
gained  in  order  to  warrant  an  additional  second  of  communication: 

profit  =  U buyer  ~  Y market  *  U seller  (7.11) 

Step  5)  Sellers  will  choose  the  bid  that  generates  the  greatest  non-negative  profit. 
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(a)  (b) 

Figure  7.5.  Experimental  market-based  task  assignment  vs.  optimal  a  priori  task  assignment 

(a)  speedup  and  (b)  time  to  completion 

In  this  way,  computational  assignments  will  be  distributed  throughout  the  network  in  such  a  way  that  the 
overall  utility  of  the  market  as  a  whole  is  maximized.  Because  of  the  addition  of  the  weighting  parameters, 
Vbuyer,  Pbuyer ,  and  Jmarket,  the  resulting  framework  is  capable  of  optimally  adapting,  in  real-time,  to 
shifting  computing  needs  within  a  wireless  network.  For  example,  assume  a  computing  task  surfaces 
where  quality  communication  channels  are  absolutely  essential.  Without  any  reprogramming  of  the 
sensing  network,  the  network  can  reassign  a  larger  Pbuyer  value  in  order  to  rcfleet  the  added  emphasis  on 
avoiding  communication  failure. 

7.3  Experimental  testbed: 

In  this  study,  four  optimization  tasks  of  varying  complexity  (the  25-Queens,  50-Queens,  75-Queens,  and 
100-Queens  problems)  are  randomly  assigned  to  four  available  Narada  wireless  sensors.  Each  of  these 
four  sensors  then  beeomes  the  “master”  node  in  the  Psa  search  tree  associated  with  their  given  //-Queens 
task.  After  these  assignments  have  been  made,  an  additional  pool  of  processing  nodes  (containing 
between  1  and  16  wireless  sensors)  is  made  readily  available  for  computational  use.  At  this  point,  the 
market-based  bidding  process  begins,  and  each  of  the  four  master  nodes  are  allowed  to  “bid”  on  the 
computational  serv  ices  of  the  free  sensing  nodes.  This  bidding  process  proceeds  as  described  in  Seetion 
7.2.3.  If  a  “master”  node  finishes  the  PSA  search  at  its  assigned  temperature  step  without  finding  a  global 
minimum,  it  w  ill  pass  its  “master”  status  on  to  its  child,  making  itself  available  for  computation  on  any  of 
the  four  computational  tasks. 

7.4  Experimental  results: 

In  order  to  evaluate  the  performance  of  the  proposed  market-based  task  distribution  methodology,  it  is 
first  neeessary  to  establish  a  benchmark  against  which  to  compare  timing  results.  In  this  ease,  in  order  for 
the  market-based  method  to  be  proven  effective,  it  must  be  shown  that  a  sensing  network  utilizing 
market-based  methods  is  eapable  of  completing  the  four  assigned  tasks  at  least  as  quickly  as  if  an  optimal 
subset  of  processors  had  been  assigned  to  each  task  at  the  outset  of  computation.  In  fact,  even  a  certain 
amount  of  degradation  in  performance  with  respect  to  an  a  priori  optimization  may  serve  to  validate  the 
method  in  this  ease,  as  the  sealable  benefits  of  real-time  task  assignment  w'ould  greatly  outweigh  a  small 
amount  of  time  savings  when  dealing  with  full-scale  problems.  For  example,  while  an  optimal  a  priori 
task  distribution  is  trivial  to  calculate  when  dealing  with  only  four  distinet  tasks,  that  optimization 
problem  grows  exponentially  as  additional  tasks  are  added. 
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Figure  7.6.  Example  record  of  market-based  task  distribution  over  time, 
compared  with  optimal  a  priori  assignment  (dotted) 

Thus,  the  first  step  in  evaluating  the  proposed  market-based  method  is  to  ealeulate  the  optimal  a  priori 
task  assignment  of  processing  units  for  each  of  the  test  eases  (5-20  Narada  units).  This  is  accomplished 
by  using  a  brutc-foree  seareh  of  the  possible  task  assignments.  Then,  experimental  data  is  gathered  using 
Narada  netw  orks  of  varying  sizes.  In  eaeh  experimental  instance,  the  network  of  sensors  is  asked  to  solve 
all  four  ^-Queens  problems.  Figure  7.5(a)  and  7.5(b)  show  the  experimental  market-based  performance 
against  the  a  priori  performance.  Figure  7.5(a)  displays  the  speedup  of  eaeh  method  relative  to  a  serial 
implementation,  and  Figure  7.5(b)  shows  the  time  required  for  eaeh  method  to  complete  the  four  tasks 
assigned  to  it.  It  ean  be  seen  from  these  plots  that  the  market-based  task  distribution  method  performs  as 
well,  if  not  better  than  an  a  priori  optimal  assignment  of  tasks.  Note  that  there  is  inherent  seatter  in  the 
market-based  results,  as  the  SA  algorithm  itself  fluctuates  widely  in  its  speed  to  convergence,  and  not 
enough  data  w'as  gathered  to  form  a  statistically  significant  eurve.  But  as  a  average,  it  ean  be  seen  that  the 
proposed  market-based  method  actually  performs  better  than  an  a  priori  optimal  distribution.  Figure  7.6 
show's  an  example  time  history  of  the  number  of  sensing  units  in  a  20-node  netw  ork  assigned  to  eaeh  task 
as  a  function  of  time.  The  dotted  lines  indicate  the  optimal  a  priori  sensor  distribution  for  a  20-node 
network  working  on  these  four  tasks. 

7.5  Section  summary  and  conclusions: 

This  section  demonstrates  a  free-market  method  of  optimally  allocating  and  consuming  scarce  system 
resources  (such  as  battery  power,  data  storage  eapaeity,  MPU  time,  wireless  bandwidth,  etc.)  within  a 
network  of  wireless  sensing  devices.  In  this  free-market  framework,  multiple  computational  tasks  are 
assigned  to  individual  sensors  w'ithin  a  w  ireless  network,  and  additional  sensors  available  for  computation 
are  distributed  amongst  these  computational  tasks  through  a  bidding  process.  In  this  seetion, 
computational  tasks  are  represented  by  simulated  annealing  seareh  problems  of  varying  complexity  Csa 
that  ean  be  broken  into  separate  bloeks  of  N$A  SA  seareh  iterations.  It  is  found  that  free-market 
distribution  methods  allow  a  set  of  multiple  computational  tasks  to  be  completed  more  quiekly  than 
possible  w  ith  any  set  distribution  of  processors,  assigned  at  the  outset  of  the  task  assignment. 
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8.  Conclusions: 

This  project  has  explored  the  deployment  of  wireless  telemetry  in  naval  plant  systems.  Unlike  traditional 
wired  counterparts,  wireless  sensors  and  wireless  actuators  can  be  installed  at  a  fraetion  of  the  cost  of 
wired  counterparts.  In  addition,  the  flexibility  of  the  wireless  eommunieation  ehanncl  supports  ad-hoc 
connectivity  whieh  allows  wireless  topologies  to  heal  when  battle  damage  oeeurs.  The  low-eost  Narada 
wireless  node  has  been  developed  eapable  of  eolleeting  data  from  scnsorS,  processing  data  at  the  node, 
aetuating  actuators,  and  communicating  with  other  nodes  via  the  wireless  ehannel. 

With  a  very  eapable  wireless  node  proposed  for  integration  in  shipboard  plants,  the  project  turned  its 
attention  to  the  computational  functionality  of  the  wireless  sensor/aetuator  network  to  perform  eondition 
assessment  (/.e.,  damage  detection)  in  the  plant.  In  particular,  model  updating  methods  were  adopted  to 
find  a  best  fit  model  of  a  plant  using  measurements  collected  from  a  network  of  wireless  sensors.  In 
particular,  simulated  annealing  is  adopted  due  to  its  computational  simplicity  and  flexibility  to  be 
distributed  in  a  network  of  computing  agents.  This  study  proposes  a  parallel  simulated  annealing 
algorithm  designed  speeifieally  to  efficiently  utilize  the  distributed  resources  available  in  large  networks 
of  wireless  deviees.  This  algorithm  gains  efficiency  as  the  number  of  sensors  in  a  network  grows,  making 
it  sealable  to  very  large  networks,  and  it  can  be  applied  to  many  of  the  large  number  of  eombinatorial 
optimization  problems  seen  across  many  engineering  disciplines.  The  proposed  algorithm  is  embedded 
within  a  network  of  Narada  nodes  to  update  an  analytical  model  describing  the  flow  of  water  in  a  chilled 
water  plant.  By  altering  the  model’s  properties  (c.g.,  pipe  diameters)  such  that  the  analytical  modal 
output  (pipe  flows)  match  experimentally  sensed  data,  properties  of  the  physical  plant  ean  be  aeeurately 
estimated.  The  method  is  proven  effective  for  estimating  rupture  conditions  in  a  chilled  water 
demonstrator  system. 

This  project  also  proposed  a  market-based  method  of  optimally  allocating  scarce  system  resources  (such 
as  battery  power,  data  storage  eapaeity,  CPU  time,  wireless  bandwidth,  etc.)  amongst  a  set  of  multiple 
computational  objectives  within  a  wireless  sensor/actuator  network.  In  this  buyer/scllcr  framework, 
available  wireless  sensors  (sellers)  are  distributed  amongst  multiple  computational  tasks  (buyers)  through 
a  utility-driven  bidding  proecss.  Because  buyers  and  sellers  in  this  market  gain  utility  in  different  ways 
(buyers  by  maximizing  speed  and  reliability  and  sellers  by  minimizing  power  consumption),  a  Pareto- 
optimal  allocation  of  scarce  resources  can  be  reached  while  completing  a  set  of  multiple  computational 
objectives  as  quiekly  as  possible.  When  evaluating  the  proposed  rcsourec  allocation  algorithm  on  a 
physical  network  of  wireless  sensor  prototypes,  it  is  found  that  this  method  allows  a  set  of  multiple 
computational  tasks  to  be  completed  as  quiekly  as  if  an  optimal  number  of  Sensors  were  assigned  a  priori 
to  each  computational  task  at  the  outset  of  computation.  This  property  is  extremely  advantageous, 
especially  as  the  number  of  computational  tasks  and/or  available  processors  increases.  By  showing  how 
this  market-based  allocation  methodology  can  be  applied  to  the  problem  of  rupture  detection  within 
shipboard  chilled  water  systems,  the  real-world  applicability  of  the  proposed  method  is  demonstrated. 


9.  Papers  Published  and  Technology  Transfer: 

A  variety  of  technology  transfer  activities  were  eondueted  in  conjunction  with  this  projeet.  In 
collaboration  with  researehers  from  the  Naval  Surfaee  Warfare  Center  Carderock  Division  (Point  of 
Contaet:  Mr.  Thomas  Brady),  the  wireless  sensors  developed  as  part  of  this  projeet  have  been  used  to 
record  the  strain  and  aeeeleration  response  of  the  FSF-1  SeaFighter,  a  high-speed  aluminum  littoral 
eombat  ship  in  current  serviee.  A  31 -ehannel  wireless  sensor  network  was  installed  to  record  hull  strains 
and  accelerations  prior  to  the  ship  departing  from  its  home  port  in  Panama  City,  Florida.  As  the  ship 
travelled  from  Panama  City  to  Portland.  Oregon,  the  response  of  the  ship  was  reeorded  to  various  sea 
states.  This  projeet  validated  the  performance  of  the  wireless  monitoring  system  on  an  operational  naval 
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vessel  while  simultaneously  showing  the  applicability  of  the  wireless  sensor  teehnology  to  hull 
monitoring  applications.  In  addition,  a  total  of  24  high  quality  publications  were  published  based  on  the 
work  conducted  in  this  project: 


1.  Zimmerman,  A.  T.,  and  Lynch,  J.  P.  “Market-Based  Frequency  Domain  Decomposition  for 
Automated  Mode  Shape  Estimation  in  Wireless  Sensor  Networks,’'  Journal  of  Structural  Control 
and  Health  Monitoring,  Wiley,  submitted [  2009. 

2.  Zimmerman,  A.  T.  and  Lyneh,  J.  P.  “A  Parallel  Simulated  Annealing  Architecture  for  Model 
Updating  in  Wireless  Sensor  Networks,”  IEEE  Sensors  Journal ,  IEEE,  in  press,  2009. 

3.  Swartz,  R.  A.  and  Lyneh  J.  P.  “Strategic  Network  Utilization  in  a  Wireless  Structural  Control 
System  for  Seismically  Exeited  Structures,”  Journal  of  Structural  Engineering ,  ASCE,  in  press , 
2009. 

4.  Lynch,  J.  P.,  Wang,  Y.,  Swartz,  R.  A.,  Lu,  K.  C.  and  Loh,  C.  H.  “Implementation  of  a  Closed- 
Loop  Structural  Control  System  using  Wireless  Sensor  Networks, “  Journal  of  Structural  Control 
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5.  Lynch,  J.  P.  and  Loh,  K.  J.  “A  Summary  Review  of  Wireless  Sensors  and  Sensor  Networks  for 
Structural  Health  Monitoring,”  Shock  and  Vibration  Digest ,  Sage  Publications,  38(2):  91-128, 
2006. 
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Wireless  Sensor  Network,”  Proceedings  of  ANCRiSST  2009:  the  5th  International  Workshop  on 
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2009. 
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Simultaneous  Execution  of  Multiple  Computational  Tasks  within  Dynamic  Wireless  Structural 
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Monitoring ,  Stanford,  CA,  2009. 
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Sensor  Network  for  Adaptive  Feedback  Control,”  Proceedings  of  the  American  Controls 
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